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1. Introduction 

The huge amount of data coming from astronomical observations in the recent years 
represents the basic impulse and the a priori condition for the development of what 
has been fairly called precision cosmology. Results on Cosmic Microwave Background 
(CMB) anisotropies by the WMAP Collaboration P and the second data release of 
the Sloan Digital Sky Survey are perhaps the most striking examples of such a 
massive effort. On one hand these experimental results give a beautiful confirmation 
of our present understanding of the evolution of the Universe and provide accurate 
information on several cosmological parameters. On the other hand they stimulated 
further efforts in increasing the level of accuracy of theoretical studies. This general 
trend is clearly recognized in the analysis of several cosmological observables, such as 
the CMB anisotropies and polarization, the Large Scale Structure formation and, for 
what matters for the present study. Big Bang NucleosjTithesis. 

Big Bang Nucleosynthesis (BBN) is one of the pillars of the cosmological model, and 
it also represents one of the most powerful tools to test fundamental physics. Presently, 
this aspect is even more relevant than in the past. In fact since the baryon density 
is now independently measured with very high precision by CMB anisotropies, the 
theory of BBN is basically parameter free in its standard formulation so that comparing 
experimental determination of light nuclides with corresponding theoretical estimate 
can severely constraint any exotic physics (for a general overview on the subject, see 
the review by B.D. Fields and S. Sarkar in the Particle Data Book [H] and references 
therein) . 

In the last decade the accuracy of BBN theory has been increased by a careful 
analysis of many of its key aspects. The accuracy of the weak reactions which enter 
the neutron/proton chemical equilibrium has been pushed up to less than percent level 
Similarly, the neutrino decoupling has been carefully studied by several authors 
by explicitly solving the corresponding kinetic equations, see e.g. [HI EI and References 
therein. These two issues are mainly affecting the prediction of ^He mass fraction, which 
presently has a very small uncertainty, of the order of 0.1 %, due to the experimental 
uncertainty on neutron lifetime. 

One of the most relevant aspects to get an accurate determination of light nuclide 
abundances is the evaluation of the several nuclear reaction rates which enter the BBN 
reaction network, as well as the corresponding uncertainties. This task involves a careful 
study of the available data or predictions on each reaction, the choice of a reasonable 
protocol to combine them in order to obtain a best estimate and an error and finally, 
the calculation of the corresponding thermal rates entering the BBN network. 

After the first analysis performed in [H1-J21! and later studies I14j. which have 
been the status of the art results for decades, an important step has been the publication 
of a large nuclear rate catalogue by NACRE Collaboration ^3] . Despite of the fact that 
this catalogue by its own nature is an all-purpose compilation and not only devoted 
to BBN studies, and that in particular not all BBN reactions are in fact covered. 
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nevertheless it is an extremely important reference. Other nuclear databases, such 
as JH] or ^Tj collect relevant information on cross sections, nuclear parameters and 
bet a- unstable nuclide lifetimes. 

Recently, new and more precise experimental data have been obtained for relevant 
nuclear reactions in the energy range of interest for the BBN, as for example the 
LUNA Collaboration results on + p ^ 7 + ^He process ^H]- In view of this new 
developments and of the fact that as stressed already there is presently a demand for an 
increased precision of BBN predictions, it is important to undertake a critical review of 
the whole BBN nuclear network, and in particular to quantify the uncertainties which 
affect the nuclide abundances as obtained by propagating the error on the several nuclear 
rates entering this network. This study is the main goal of this paper. Similar analysis 
has been also performed in ^HI, jSOj and |21j . 

In the following Sections we cover all the main theoretical aspects of BBN with 
the aim of providing a self contained description of the several issues which enter an 
accurate determination of nuclide abundances. To help the reader to distinguish the 
review parts collecting the results of previous studies from the ones devoted to original 
new results first presented in the present paper, the former have been tagged by a "*" . 

In Section 2 we review in details the general BBN framework by recalling the 
corresponding set of equations. We also describe the neutrino decoupling as obtained 
via a numerical solution of the corresponding kinetic equation, including as in f?! the 
effect of QED radiative corrections and in particular the non thermal distortion in 
neutrino distribution functions due to e"*" — e~ entropy release. 

Section 3 is devoted to the analysis of all main rates entering the BBN network. We 
first review the calculation of <-> p weak rates, evaluated up to order a QED radiative 
corrections, and also taking into account the effect of finite nucleon mass and thermal 
radiative effects. We also discuss the effect on these rates of neutrino decoupling. After 
introducing the nuclear astrophysics formalism used in our analysis, we describe our 
method for the data reduction and rate estimate and finally describe in details the 
data and the results obtained for all leading reactions entering the BBN network. We 
also analyze some sub-leading reaction which, in view of the present uncertainties on 
the corresponding rates, may still play a role in determining the eventual light nuclide 
yields. 

In Section 4 we first summarize the present status of observations on ^H, ^He, and 
''Li, as well as what is known on other light nuclides, such as ^He or ^Li. We then give 
our results for main nuclide abundances, as obtained by a FORTRAN code that we 
have developed in the recent years modifying the original public code of [22] in order 
to consistently introduce all the issues described in this paper (neutrino decoupling, 
radiative corrections ion ^ p rates, etc.) and employing a different numerical resolution 
method, as explained in Section 2.3. 

All results are given for a standard BBN scenario, with three non degenerate 
neutrinos, and the baryon density as fixed by the WMAP result combined with results 
of CBI and ACBAR experiments on CMB anisotropies, as well as with 2dFGRS data 
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on power spectrum, cj;, = = 0.023±0.001 PP, where is the ratio of the baryonic 
matter density with respect to the critical one, and h is the Hubble constant in units 
of 100 Km Mpc~^. A careful study of this case seems to us of particular relevance 
since there are no free parameters left. Comparing experimental results with theoretical 
expectations is therefore a clean consistency test of the simplest BBN dynamics. It may 
also give useful hints to point out possible systematics in the experimental results for 
nuclide abundances. A particular stress is given to the present theoretical uncertainty 
on each nuclide yield, as well as on the role of the several nuclear rates in building up 
this uncertainty. We also comment on the impact of different values adopted for uoi, as 
obtained by WMAP with different choices of priors or combining different observations. 
Finally, we discuss how the baryon density is constrained by BBN as compared to the 
very precise measurement by CMB anisotropics. 
Our Conclusions are reported in Section 5. 

2. The BBN general framework * 

2.1. Generalities * 

We consider Nnuc species of nuclides, whose number densities, n^, are normalized with 
respect to the total number density of baryons, n^, 

X, = — 2 = n, p, ^H... . (2.1) 

riB 

A list of all nuclides which are typically included in a BBN analysis is reported in Table 
1. To quantify ^H, ^He and ^Li abundances, we also use in the following the parameters 

X2h/Xp, XsHe/Xp, Yp = AX4iie, Xr^/ Xp , (2.2) 

i.e. the ^H, ^He and ^Li number density normalized to hydrogen, and the ^He mass 
fraction||. For the sake of clarity we stress that our notation Xi for the nuclide number 
density normalized to hydrogen is not adopted by all authors. The symbols Yi can be 
also frequently found in the literature. 

The neutrino/antineutrino distributions are denoted by fu^dpl ,t), /p^dpl 7^) and 

U=U^UMp\,t) , f,^ = f,^^f,^(\p\,t) , (2.3) 

In the temperature range we are interested in, lOMeV > T > 0.01 MeV, electrons and 
positrons are kept in thermodynamical equilibrium with photons by fast electromagnetic 
interactions. Thus, they are distributed according to a Fermi-Dirac function /e±, with 
chemical potential fie- 

II Though this definition is widely used, and we will use it as well, yet it is only approximately related to 
the real mass fraction, since the ^He mass is not given by 4 times the atomic mass unit. The difference 
is quite small, of the order of 0.5%. However, in view of the present precision of theoretical analysis on 
^He yield, this difference cannot be neglected and we think it is worth to stress this point in order to 
avoid possible misinterpretations. 
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1) n 


7) ^Li 


13) 


19) 


25) ^^0 


2) p 


8) ^Li 


14) 


20) 


26) 


3) 2H 


9) ^Be 


15) 


21) 




4) 


10) ^Li 


16) i^B 


22) 




5) =^He 


11) SB 


17) 


23) 




6) ^He 


12) ^Be 


18) 


24) 





Table 1. Nuclides considered in the BBN analysis 



The set of differential equations ruling primordial nucleosynthesis is the following 
(see for example 



— = -3i/ , (2.5) 

riB 

p = -3H{p + p) , (2.6) 




(2.8) 

/ ^- — 

with Pa = ^e, t'e, i^x, i^x, where p and p denote the total energy density and pressure, 
respectively, 

p = py + pe + P,. + PB = PNB + PB , (2.10) 
P = P-f+Pe + Pu + VB = PNB+VB , (2.11) 

where i,j,k,l denote nuchdes, pb and pnb are the baryon and non baryonic energy 
density respectively, Zj is the charge number of the i—th nuclide, and the function 
L{C,,y) is defined as 

dCC\/C^(j^-^) . (2.12) 

Equation ()2.4|1 is the definition of the Hubble parameter, H, with Gn the gravitational 
constant, whereas Equations ()2.5|) and ()2.6|) state the total baryon number and entropy 
conservation in the comoving volume, respectively. The set of Nnuc Boltzmann equations 
fl2.7|) describes the density evolution of each nuclide specie, with T^i^ij the rate per 
incoming particles averaged over kinetic equilibrium distribution functions. While in 
fact chemical equilibrium among nuclides cannot be assumed, as BBN strongly violates 
Nuclear Statistical Equilibrium, it is perfectly justified to assume kinetic equilibrium, 
which is maintained by fast strong and electromagnetic processes. Equation ()2.8j) 
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states the Universe charge neutrahty in terms of the electron chemical potential, with 
= ^e/T and T the temperature of e^,7 plasma, and finally Equations ()2.9|) are the 
Boltzmann equations for neutrino species, with 1^^ [fu^, fyj\ standing for the collisional 
integral which contains all microscopic processes creating or destroying the specie u^- 
In the following we do not consider neutrino oscillations, whose effect has been proved 
to be sub-leading and affect ^He mass fraction for ~ I-IO""^ |26j . 

In the standard scenario of no extra relativistic degrees of freedom at BBN epoch 
apart from photons and neutrinos, the neutrino chemical potential is bound to be a 
small fraction of neutrino temperature, smaller than approximately |^,^| = \^y/Ty\ < 0.1 
|2() | n7 \ I29j. This bound applies to all neutrino flavors, whose distribution functions 
are homogenized via flavor oscillations [2ZlEni- In view of this, as far as this analysis is 
concerned, we will focus on non degenerate neutrinos, so that fi,^ = fp^ and fu^ = fp^. 

The neutrino energy density and pressure are defined in terms of their distributions 

as 

p. = 3 p. = 2 I ^ \p\ + 2 U , (2.13) 
whereas for baryons we have 



PB 



i 

with the AMi and M„ the i-th nuclide mass excess and the atomic mass unit 
respectively. Finally, pressure and energy density for the electromagnetic plasma (e^ 
and 7) is calculated by improving the incoherent scattering limit result (particles 
correspond to poles of the propagators at, respectively, = (photons) and p"^ = ml 
(electrons/positrons)) and considering the effect of finite temperature QED corrections. 
It has been shown in fact that at the time of BBN, these corrections affect Equations 
fl2.4|) - ()2.9|) at some extent and slightly influence the '^He abundance j3]. Indeed finite 
temperature QED corrections modify the electromagnetic plasma equation of state and 
thus influence Equation ()2.6p . the expression of the expansion rate H and, as it will be 
discussed in the following, the conversion rates n ^ p. 

The change in the electromagnetic plasma equation of state can be evaluated by 
considering the corrections induced on the and photon masses, using the tools of 
Real Time Finite Temperature Field Theory For the electron/positron mass, up 
to order a = e^/ (4 7r), we find the additional finite temperature contribution [71021 

5ml{\plT) = + — / d\k\ 

6 TC Jo 



2 mi a „, , Ik 



/ ^ikiyiog 

Jo Ek 







1 






IpI - 


f |k| 


1 


IpI 


- |k| 





(2.15) 



where = y |kp + m"^. We note that in this equation the term depending on the 
momentum |p| contributes for less than 10 % to Sm^ and thus it can be safely neglected. 
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This highly simphfies the analysis and is fully legitimate at the level of precision we are 
interested in (corrections to Yp at the level of 0.0001). 

The renormalized photon mass in the electromagnetic plasma is instead, up to order 
a, given by (see for example ^33i) 

Sm' (T) = — ^|k|^ ^ . (2.16) 

^ n Jo Ek e^fc/^ + 1 

The corrections (|2.15j) and (j2.16p modify the corresponding dispersion relations as 



|k|2 + ml + 5mj (T) 



e, 7). Thus the total pressure and energy density 



of the electromagnetic plasma result to be 



J" roo 

P7 + Pe = ^ / d\k\ |k|^ log 

TT^ JO 



1 + 

^1 _ Q-E,/T^ 



P-l + Pe 



T 



dT 



1 (Po 



Pe 



(2.17) 



(2.18) 



Expanding these expressions with respect to 5ml and 5m^, one obtains the first order 
correction 

|kM|k| 



5p^ + (5pe 



2 7r2 



|k| 5ml {T) 1 (5m^ (T) 
1 + 2 e^/T - 1 



(2.19) 



Eu e^'^/^ 

The energy density is then obtained by using the expression for p^ + pe in Equation 
(ITTHll . 

The equations ()2.4p . ()2.6|) and ()2.9|) can be solved separately since they are 
essentially determined by relativistic particles, and only negligibly affected by the 
baryons contribution^. This allows to solve the evolution of the neutrino species first, 
and then to substitute the result into the remaining equations. 



2.2. Neutrino decoupling * 

The evolution of neutrino distribution function can be highly simplified by using the 
scale factor as evolution variable x = mg a, and the comoving momentum y = \p\ a. We 
also define the rescaled photon temperature as ^ = T a. By using these definitions, the 
set of equations ()2.6|) and ()2.9|) . once baryonic contribution is neglected, becomes 



f'^ci (-^J y) ~ FT '^'^a [fi^e 5 ffx 



with z/„ 



(2.20) 
(2.21) 



dx """"^ J ^ / 

In Equation (|2.20p . which states the conservation of the total energy momentum, p and p 
are the dimensionless energy density and pressure of the primordial plasma, respectively. 



X 

P = P{ — 

mp 



X 

Pnb I — 

me 



P = p(— )'~Pjvb(— . (2.22) 



^ The effect of baryons on neutrino decoupling is due to nucleon-neutrino scattering processes, as well 
as to the baryon contribution to Hubble parameter pB- Both effects are expected to produce very small 
corrections of order 77 or M^^rj/T. 
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As in Reference [71 1^ the unknown neutrino distributions are parameterized as 

1 



fv^ (a;, y) 



1=0 



where Pi{y) are orthonormal polynomials with respect to the Fermi function weight 
dy 



(2.23) 



/„ e. + i • (^■^*) 

With these definitions the problem of finding the (momentum dependent) distortion 
in neutrino distribution function is then reduced to determine the time evolution of 
suitable linear combinations of the lower momenta of these distributions. This method 
has been shown to provide results in quite a good agreement with different approaches, 
which instead use discretized comoving variables and solve Boltzmann equations on a 
grid in the y variable [SSI IHEj • 

By substituting Equation ()2.23|) into Equations ()2.2H) . including the QED 
corrections discussed previously, and using the covariant conservation of the energy 
momentum tensor (|2.2(jp as an evolution equation for one gets 

d-z _ i-A^I^)- A C dyy^ + 2^) + G,(x/. 
dx 
d 



^J{x/z) + Y{x/z) + ^ + G2{x/z) 



where 

Gi{uj) = 2na 
'K'(u) 



dx ^ ^ ' xH 
"1 (K{u 







dyPiiy) [U^Ju^ 



(2.25) 
(2.26) 



J(^) 



UJ 



3 +^*'<-)'- 6 



- K{lj) J{uo) 



+ 



6 



Kiuj) K'(uj) + 



J'(co) 



6 



G2{uj) = -Sna \ 



'K{uj) J{uj) 



6 



6 



J'{uj) K{uj) + J{uj) K'{uj) 



-K{uof + K{uo) J{uo) 



(2.27) 



with 



6 



- K{uj) K\uj) + 



+ J\uj) K{uj) + J{uj) K\iu) 



TT^ JO 



du 



U 



yju^ + uj"^ exp [y/u^ + ujA + 1 



J{u) 



TT^ Jo 



du u 



du u 



exp 



(exp (Vu2 + Cj2) + 1 
exp 



(2.28) 

(2.29) 
(2.30) 
(2.31) 



(exp (Vu2 + + l) 

The functions K'{uj) and J'{oj) stand for the first derivative of K{(jj) and t/(co') with 
respect to their argument. 
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Figure 1. The evolution of z 
small temperatures is z= 1.3984 



T/Tn versus z — ni^/T. The asymptotic value at 



Since at high temperature, T ~ 10 MeV, the neutrinos are in thermal equihbrium 
with the electromagnetic plasma, the initial condition for the coefficients is a" = and 
we can always choose the initial value for the scale factor such that we have, at 10 MeV, 
z = 1. From Equation ()2.25|) . neglecting the terms proportional to the derivative 
of neutrino distributions, as well as QED corrections, one gets the asymptotic value 
z^q = (11/4)^/^ = 1.4010 which represents the ratio between the photon and neutrino 
temperatures after the complete annihilation of e~^e~ pairs into photons only. This limit 
is often referred to as the instantaneous decoupling value. A fully numerical treatment 
of the decoupling shows indeed that the asymptotic value of T/T^, z = 1.3984, is 
lower than since neutrino plasma is slightly heated by the e~^e~ annihilations. We 
report in Figure [T] the actual behavior of z. In [71 E3j it was shown that a very good 
approximation, giving results with accuracy of 1%, is to consider polynomials in ()2.23p 
up to third order. The neutrino distribution function can be also written as 



In Table 2 we report the asymptotic values of the c^ coefficients, as found in 17], while 
their evolution versus z is shown in Figures |21 and 121 respectively. 

2.3. Numerical solution of the BBN set of equations * 

Once neutrino distribution functions are determined, the BBN set of equations is reduced 
to Equations (2.4)-(2.8) of Section 2.1. This set can be further reduced. Actually it is 
more convenient to follow the evolution of the Nnuc + 1 unknown functions (0e, Xj) in 




(2.32) 
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Figure 2. The evolution of the electron neutrino distortion coefficients Cq (solid), cf 
(dashed), cl (dotted) and cf (dot-dashed) versus z = nie/T 
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Figure 3. The evolution of the fi,T neutrino distortion coefficients eg (solid), cf 
(dashed), cf (dotted) and cf (dot-dashed) versus z = irie/T 
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Flavor (a) 




cf 






e 


-2.507 


-2.731 


6.010 


-0.1419 


/i,r 


-2.003 


-2.196 


3.061 


-0.1091 



Table 2. Values of the coefBcients of Equation H2.32|l in unit of 10 ^ 



terms of the dimensionless variable z = rUe/T, and to use Equation ()2.8|) to get ub as 
a function of 0e- The new set of differential equations may be cast in the form 
#e _ 1 L E F + {zL,-3L) G 
dz z 



L E ^ 



dpe 



Lrh^. G 



(2.33) 



dX, 
dz 



z 



L^^F + {zL,-3L) II 
L E ^ - L^G 



(2.34) 



where the functions E, F and G are given by 
E{z, (l)e,X,)=3H-^ 



J2 ^j^j 



F{z 



dpe dp^ 



— z- 



' 2 ~ dz dz 



G{Z, 0e, Xj) =3H [ Pe,^ + Pe,^ + Pb + 



Niz) 



+ 



zL 



J2j Zj Xj 



E AM. 



2z 



r, 



with 



N{z) 



z^ 



' d ^ 
. dx 



(2.35) 
(2.36) 

(2.37) 
(2.38) 



x=x{z) 

and H = nieH, ub = rn^ fiB, Fj = m^Ti, p = T'^ p, p = T^p, and finally 
Pu = Puix/meY. Note that by using the previous definitions, it is possible to express 
Pb and p^ as functions of z, (p^, and Xj only 



z Liz, 



Pb 



Pb 



J^j Zj Xj 



m. + E(am, + — j X, 



(2.39) 
(2.40) 



With AMj and M„ we denote the i-th nuclide mass excess and the atomic mass unit, 
respectively, normalized to mg. The values of the partial derivative of L with respect 
to z and 0e, denoted with and L,^^, and the quantities pe, dp^/dz and (9pe/90e in a 
form which is suitable for a BBN code implementation can be find in Appendix A of 
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The neutrino contribution to the previous equations, via the function N{z), can be 
obtained by the solution of Equations ()2.25|) . ()2.26|) . written as function of the variable 
z. To this end it is necessary to invert (numerically) the relation 

= z — >x = x{z) . (2.41) 

z{x) 

It is interesting to notice that neutrinos only contribute to the z evolution of 0e and Xj 
via the non equilibrium terms in their distribution functions, since 

d A / 1 /■ , 9^,- 1 \ d 



X 



This expression can be numerically calculated as a function of x, and then it should be 
evaluated at x{z). This result is quite expected. In fact neutrinos well before decoupling 
share the same temperature of the electromagnetic plasma, so z is a constant in this 
limit. On the other hand, once they are fully decoupled at low temperatures, they satisfy 
an independent entropy conservation condition, and so they do not further contribute 
to the time evolution of z. Only during the electron/positron annihilation phase they 
do affect this evolution, via the non equilibrium terms 5fi^^{x,y), which in fact are the 
genuine effects of their residual coupling to the electromagnetic plasma. 

Equations ()2.33|) - ()2.34|) are solved by imposing the following initial conditions at 
/(lOMeV) 

M^in) =0e° , (2.43) 
XniZin) = {exp{q Zin} + , (2.44) 

Xp{zin) = {exp{-q Zin} + l)'^ , (2.45) 
^(...)^|(c(3,/f)^'"^|(£-)''^-%-X^X- 

X exp {BiZin} I = 'H, ^H,... . (2.46) 

In the previous equations q = {Mn — Mp)/me, and the quantities Ai and Bi denote the 
atomic number and the binding energy of the i— th nuclide normalized to electron mass, 
respectively. Finally rji is, the initial value of the baryon to photon number density ratio 
at T = 10 MeV (we will discuss the final to initial value ratio of rj in the following), and 
the solution of the implicit equation 

L{z,n, 0e°) = V^ E Mz^n) ■ (2.47) 

i 

The method of resolution of the BBN equations ()2.33p . ()2.34p is the same apphed 
in [21]. It uses a class of Backward Differentiation Formulas with Newton's method, 
implemented in a NAG routine with adaptive step-size (see |24i for more details). As 
an example of the evolution of nuclei abundance we report in Figure ID the result of our 
numerical code for Uh = 0.023 and three neutrinos versus m^/T. 
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Figure 4. Nuclear abundances depart from their equilibrium values, undergo an out- 
of-equilibrium phase, and eventually reach their final values 



3. The Nuclear rates 
3.1. Weak Reactions * 

The weak reactions transforming n ^ p, namely 

(a) z/g + n — > + p , (d) Ve + p ^ + n , 

(b) e" + p Ue + n , (e) n ^ e~ +T'e + P , 

(c) e'^+n^Ue+p , (/) e'+V^+p—^n . (3.1) 

are the leading processes in fixing the neutron abundance at the onset of BBN and thus 
a key quantity in determining the ^He mass fraction. In view of this, much effort has 
been devoted to refining the theoretical accuracy in evaluating these processes, which 
presently is at the order of 0.1%. We here summarize the main results obtained in the 
literature, while more details can be found in IH| I24j . 

The Born rates are the tree level estimates obtained with V — A theory and with 
infinite nucleon mass. As an example, for the neutron decay process (e) we have 

ujBin^e +Ve+p) = d\p'\\p\'q^e{qo)x 

X [l-f,^{q,)][l-f^{p'^)] , (3.2) 

where Gp is the Fermi coupling constant, Cy and Ga the nucleon vector and axial 
coupling. According to our notation p' and p'q are the electron momentum and energy, 
and go = — Mp — p'^ = A — p'q the neutrino energy. The rates for all other processes 
(a) — (d), (/) can be simply obtained from (|3.2|) properly changing the statistical factors 
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and the expression for go (sse for example [3]). Average is performed at this level of 
approximation over equilibrium Fermi-Dirac distribution for leptons, i.e. neglecting the 
effects of distortion in neutrino/antineutrino distribution functions. In Figure El we 
report the Born rates ujb for n ^ p processes, The accuracy of Born approximation 




-7- 

10 r 



T/m. 



Figure 5. The total Born rates, ub, for n — > p (solid line) and p n transitions 
(dashed line). 

results to be, at best, of the order of 7%. This can be estimated by comparing the 
prediction of Equation ()3.2|) for the neutron hfetime at very low temperatures, with the 
experimental value r^^' = (885.7±0.8) s 

A sensible improvement is obtained by considering four classes of effects 

1) electromagnetic radiative corrections, which largely contribute to the rates of the 
fundamental processes; 

2) finite nucleon mass corrections, which are of the order of T/Ma? or rrie/MTv, M^v 
being the nucleon mass; 

3) thermal/radiative effects, proportional to the surrounding plasma temperature; 

4) non instantaneous neutrino decoupling effects on the Born rates. Distortion in 
electron neutrino distribution function directly enters the thermal average of weak 
rates, as well as the detailed evolution of the temperature ratio T /Ty = z. 

3.1.1. Electromagnetic radiative corrections. Electromagnetic radiative corrections to 
the Born amplitudes for processes ()3.1|) are typically split into outer and inner terms 
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(for a review see e.g. [SZl)- The first ones involve tlie nucleon as a whole and consist of 
a multiplicative factor to the squared modulus of transition amplitude of the form 

1 + :t-^?(Po,9o) . (3.3) 
/vr 

The function g{pQ,qo), can be found in Reference jHHl, and depends on both electron 
and neutrino energies. On the other hand, the inner corrections are deeply related to 
the nucleon structure. They have been estimated in Reference j^HI; and result in the 
additional multiplicative factor 

where the first term is the short-distance contribution and Ag = —0.34 is a perturbative 
QCD correction. The other two terms are related to the axial-induced contributions, 
with Ma = (400 ^ 1600)MeV a low energy cut-off in the short- distance part of the 'jW 
box diagram, and C is related to the remaining long distance term. 

Summing up these two kind of corrections, and resumming all leading logarithmic 
corrections a"/n"(M^) ^j, one gets the global multiplicative factor 




27: \ Ma J 1-n 



S{Mp,Mz) , (3.5) 



where a(/i) is the QED running coupling constant defined in the MS scheme and 
S{Mp,Mz) a short distance rescaling factor, defined in [S]. 

Another electromagnetic effect is present when both electron and proton are either 
in the initial or final states (namely, processes (a), (6), (e) and (/)). It is the Coulomb 
correction due to rescattering of the electron in the field of the proton and leading to 
the Fermi function 



Tip',) ^ 1 + . (3.6) 




Summing all these corrections and using the most recent estimates for G^? = 
(1.16637±0.00001)-10"^GeV"^ Cy = 0.9725±0.0013 and the ratio Ca/Cv = 
— 1.2720±0.0018, |3|, and including the (small) finite nucleon mass corrections discussed 
in the next Section, the theoretical prediction for neutron lifetime, which as we said can 
be used as a check of accuracy, is now = 886.5 s. The inclusion of electromagnetic 
radiative corrections therefore gives quite a satisfactory result for r^, with accuracy of the 
order of 0.1%. As in |211E1 we deal with this tiny residual discrepancy by assuming the 
presence of a further multiplicative factor which properly renormalizes the theoretical 
prediction. Such a factor 1 + 6r = t^^/t^^ = 1.001, is then applied as a rescaling factor 
to all processes ()3.1|) . Notice how this value is remarkably closer to unity with respect 
to what estimated in [2ll IH] . This is simply because of the different experimental value 
of r^^ quoted in |3] with respect to previous estimates. In closing we have to mention 
that there may be still a systematic effect induced by adopting this procedure, since it 
assumes that all residual discrepancy parameterized by Sr is independent on energies 
of the outgoing particles. However it is worth saying that with present data on "^He 
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abundance, this effect is hardly playing a role in a comparison of theoretical prediction 
for BBN with data. 

3.1.2. Finite nucleon mass corrections.* There are relevant contributions to the n ^ p 
rates which appear by relaxing the infinite nucleon mass approximation. These effects 
are proportional to mf.lMfq or T /M^q, and in the temperature range relevant for BBN, 
can be as large as the radiative corrections. This has been first pointed out in j3T] and 
then also numerically evaluated in |U E] . 

At order 1/Mjv, the weak hadronic current receives a contribution from the weak 
magnetic moment coupling 



where, from CVC, /2 = — /^n)/2 = l.SlVud- Both scalar and pseudoscalar 

contributions can be shown to be much smaller and negligible for the accuracy we are 
interested in. At the same order l/M^ the allowed phase space for the relevant scattering 
and decay processes get changed, due to nucleon recoil. Finally, one also has to consider 
the additional contribution due to the initial nucleon thermal distribution. In fact in 
the infinite nucleon mass limit the average of weak rates over nucleon distribution is 
trivial, since the nucleon is at rest in any frame. For finite M^r, by considering only 
1/Mn terms, the effect of the thermal average over the thermal spreading of the nucleon 
velocity produces a purely kinetic correction Aujk, whose expression can be reduced to 
a one-dimensional integral over electron momentum and then numerically evaluated. 
The explicit expression, which we do not report for brevity, can be found in Section 4.2 
and Appendix C of 0. 

3.1.3. Thermal- Radiative corrections. The n ^ p rates get slight corrections from the 
presence of the surrounding electromagnetic plasma. To compute these corrections once 
again one may use the standard Real Time formalism for Finite Temperature Field 
Theory to evaluate the finite temperature contribution of the graphs reported in 
Figure El for the n p processes. Inverse processes p ^ n are obtained by inverting 
the momentum flow in the hadronic line. The first order in a is given by interference 
of one- loop amplitudes of Figure El b) and c) with the Born result (Figure El a)) . As 
usual, photon emission and absorption processes (Figure Eld)), which also give an order 
a correction, should be included to cancel infrared divergences. Notice that photon 
emission (absorption) amplitudes by the proton line are suppressed as . 

All field propagators get additional on shell contributions proportional to the 
number density of that particular specie in the surrounding medium. For 7 and e^, 
neglecting in this case the small electron chemical potential, we have 




(3.7) 




(3.8) 




(3.9) 
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Figure 6. The tree level Born (a), the one-loop (b),(c), and the photon 
emission/absorption diagrams (d) for n ^ p processes. 



with the photon distribution function. The entire set of thermal/radiative corrections, 
at first order in its typical scale factor, i.e. aT/m-e, have been computed by several 
authors jl2]-jnij with quite different results. 

It was correctly pointed out in [55^ that when order a QED corrections are 
introduced there are new processes taking place in the plasma which should be added 
and contribute to neutron/proton chemical equilibrium 

7 + p— ^e^ + z/g + n , e'^ + i'e + n^'^ + p . (3.10) 

For completeness, since the processes flH.lOj) where not considered in we report the 
explicit expressions for their rates, evaluated in the infinite nucleon mass limit 

Gl {Cl + ?,Cl)a 



^(7 + p e'^ + i>e + n) 



47]-3 



TT JO 



Ik'l Ikl 



X 



/co (A + fco + |k'| 



d|k' 

(A + A;o + |k'|)' 



ci|k| 



k'l 



log 



^0 + 


|k| 


ko - 


|k| 



(3.11) 



+ 2 |k'| (A + |k'|)] /,(A + ko+ |k'|) (1 - /,+ (A;o)) (1 - ^(|k'|)) 
where ko = a/ |kp + m^. Inverse process is simply obtained by changing the statistical 



(3.12) 



factor with the following one 

(l + /,(A + A:o + |k'|))/e+(M/.e(|k1) . 

The relative corrections to the total rates are shown in Figure U\ 

In the range of temperatures when neutron fraction freezes out, this rate is not 
severely suppressed by conservation of energy-momentum, since photon mean energy is 
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of the order of MeV. Nevertheless it should be pointed out that the freeze out of neutron 
to proton ratio is mainly dictated by two body processes (a) — {d) , which in fact dominate 
over neutron decay and inverse process (e) and (/) in the relevant temperature range, 
so the inclusion of ()3.10|) is very weakly affecting Yp. In the following we will adopt 
the results for thermal corrections obtained in Reference [Sj, to which we refer for all 
details, updated with the inclusion of the ()3.10p and inverse process contributions. The 




Figure 7. The relative correction to the n ^ p (solid line) and p ^ n (dashed line) 
total rates due to the ^+p^e^ + h'e + n and inverse processes. 

effects of the QED interaction in the plasma also show up in modifying the equation of 
state of electron/positron and photons, as well as they induce thermal contribution to 
their mass. These effects, first considered in jj], have been already described in Section 
12.21 In particular the thermal averaged weak rates get an order a corrections when the 
electron thermal mass is introduced in the corresponding distribution function. 

3.1.4- Non instantaneous neutrino decoupling effects on the Born rates. As we have 
already described, neutrino distribution functions get distorted for the effect of a 
partial entropy release to neutrinos by electron/positron pairs. Actually, the e"*" — e~ 
annihilation phase is not an instantaneous phenomenon, but it partially overlap in time 
with n/p ratio freezing, as discussed in jSH]. The overall effect on BBN is due to three 
different phenomena. First, the weak rates (a) — (/) are enhanced by the larger mean 
energy of electron neutrinos. On the other hand there is an opposite effect due to the 
change in the electron/positron temperature. Finally, since the photon temperature is 
reduced with respect to the instantaneous decoupling value = (11/4)^/^, the onset of 
BBN, via synthesis, is taking place earlier in time. This means that fewer neutrons 
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decay from the time of freezing out of weak interactions and this in turn corresponds to 
a larger ^He yield. 

All these effects have been considered in our analysis. The actual behavior of 
shown in Figure ^ has been used in the equation ruling BBN, as described in Section 
12. HI Moreover the Born rates for weak processes are corrected by using the distortion 
^ f 1 y) 1 computed numerically. The resulting behaviors of the relative corrections 
to the total n ^ p rates is shown in Figure |H1 In particular 6uj is the difference of 
the n ^ p rates as from present analysis and the result obtained in the instantaneous 
neutrino decoupling limit, with decoupling occurring at = 2.3 MeV. Notice the 
discontinuity in the first derivative at Tq. In fact, d(Ty/T)/dT has a discontinuity 
in the instantaneous decoupling limit, while it has a smooth behavior when neutrino 
decoupling is studied via a full numerical solution of the corresponding kinetic equations. 
We find that the effect of distortion of neutrino distribution function slightly increases 
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Figure 8. The relative correction to tlie n ^ p (solid line) and p ^ n (dashed 
line) total rates, due to the neutrino distortion. Suj is the difference of the n p 
rates calculated by numerically solving the neutrino decoupling and the ones for an 
instantaneous neutrino decoupling at To = 2.3 MeV. For p ^ n the effect is shown 
down to T ~ 0.3 MeV, since for lower temperatures the total rate becomes negligibly 
small. 

the Yp, by the small amount Al^ ~ 1-10~^, in agreement with the result of |56j . 

3.1.5. The total rates for n ^ p reactions. Apart from the effect of neutrino distortion 
and of the p + 7 process discussed in the previous Sections, the effect of all remaining 
corrections listed in Section 3.1 to the weak rates has been discussed in details in jUinj. 
We here report the main results, for the sake of completeness. The leading contribution 



n ^ p 

- p ^ n 




± 
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is given by electromagnetic radiative corrections, which decrease monotonically with 
increasing temperature for both p ^ n and n —>■ p processes. Including also the effect 
of finite mass corrections due to weak magnetism and 1/Mn corrections to phase space 
the effect ranges from -2 % to 8 % for n ^ p and -3 % to 7 % for p — >• n, respectively. In 
particular they affect the rates for a 2% at T ~ IMeV, i.e. at the freeze out temperature. 
Finite mass corrections Aujk (see Section 3.1.2) are less important and contribute for 
at most 0.5 % ^ 1.5% again in the temperature range where the freezing phenomenon 
takes place. Finally the effects due to plasma corrections and thermal radiative effects 
are sub-leading, changing the rates at the level of (0.3 -r- 0.6)% only. 

In order to use the n ^ p rates in a numerical code, it is useful to fit their expressions 
as a function of z. The result is reported in [Appendix C The fit has been obtained 
requiring that the fitting functions differ by less than 0.01% from the numerical values. 
Notice that it is also a good approximation to consider a vanishing rate a;(p — >■ n) for 
T < 0.1 MeV, see Equation ()C.2|) . since it is a rapidly decreasing function with T — > 0. 

As a check of the result it is interesting to look whether the n ^ p rates we have 
obtained, satisfy the detailed balance condition, namely 

uj(n v)n„ 

p{T) = ^ ^ = 1 , 3.13 

where n„(p) are the neutron (proton) number density. In Figure 5 we plot the ratio 
p(T) versus temperature in unit of rUe- We also show the behavior of z^^ . The result 
shows indeed that detailed balance condition is satisfied with very good accuracy for 
T > me, better than the 1% level, smaller than the typical order of magnitude of order 
a radiative corrections in the same temperature range, which set the level of accuracy 
of our results. The steep rising of p{T) for smaller temperatures, T < me is easily 
understood, since for T < rUe pair annihilation phase starts and neutrino do not share 
anymore the same temperature of the electromagnetic plasma, as can be seen by the 
behavior of z~^ = T^/T. Neutrinos become colder than so that all processes with 
an ingoing neutrino/antineutrino are suppressed with respect to the reverse ones. In 
particular, the detailed balance condition for the neutron decay, the leading neutron 
destruction channel for temperatures lower than 0.3 me, and the inverse e~ + i^e+P ^ ^ 
reaction does not hold anymore and this fact leads to the increase of p(T). Stated 
differently, this large departure at low temperatures of p{T) is simply the effect of the 
lack of thermodynamical equilibrium due to the freezing of weak interactions. This 
temperature region is however irrelevant for neutron to proton ratio which has been 
already frozen out well before the annihilation phase. 



3.2. Nuclear reaction rates: analytical calculations. 

In the primordial plasma during BBN the energy of colliding nuclei is always much 
lower than their masses, and their density is very far from values needed to produce 
relevant quantum effects. It follows that the averaged reaction rates (av) for a typical 
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Figure 9. The detailed balance ratio p{T) (solid line) versus temperature in unit of 
rUe- The dashed line is the neutrino to photon temperature ratio z~^. Dotted line 
corresponds to p{T) for the case = T. 



a + h reaction process can be obtained by folding the cross section cr(-E') with a Maxwell- 
Boltzmann phase space distribution 

i^av) = . r dEEa{E) e-^/^ . (3.14) 

with fiab the reduced mass of a and b and E denoting the kinetic energy in the Center of 
Mass (CM) frame. In the following Sections we will perform this integral for the cases 
of interest for the BBN. Though this issue can be found in many excellent manuals, 
jnZl-|nni; we wm briefly summarize the main aspects of it, for the sake of fixing notation 
and to summarize the parameterizations adopted for reaction rates, which in some case 
have been improved with respect to standard results. 

3.2.1. Non-resonant reactions induced by neutrons.* Assuming a s-wave collision, as 
it is typically the case for low energies and far from resonances, for a neutron-induced 
reaction one finds a oc v~^, from which av ~ constant. More generally one defines the 
auxiliary function 

R{E) ^ icrv){E) ^ a = ^ = ^R{E) , (3.15) 
in terms of which, the rate can be written as 

M = y-y^ dy^R{yT)e'^ , (3.16) 
where y = E/T. 
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For a non-resonant reaction i? is a weakly dependent function of energy E, which 
allows to approximate R{E) by a low order polynomial in E^^'^ 



R^E) = J2^^^E^^' , (3-17) 

n=0 ^■ 



getting 



where 



{av){T) = J2bnT-/' , (3.18) 

n=0 

3.2.2. Non-resonant reactions induced by charged particles. As well known, the 
astrophysical S factor is defined as 



SiE) = aiE)Eexp[^EG/E] , (3.20) 

with Eg = 2'}T'^fiab{ZaZhay the Gamow energy. As in the case of neutrons, for a non- 
resonant cross section S{E) is a slowly changing function, and thus can be typically 
parameterized as a (low-order) polynomial approximation in E. In terms of S{E), one 
can write 



{av) = J^T-^/^ r dyS{yT)e-'^^y'^^^ = J^T-^/^I , (3.21) 

V nflab Jo V ^l^ab 

where y = E/T, yG = Eg/T and (f){y, yG) = y + Expanding (f){y, yG) up to second 

order around its minimum at Gamow peak, y^ = {yc/^^V^^ and assuming that, in the 
neighborhood of yo, S{E = yT) c::^ S{Eq = y^T) = Sq, one gets 



{av)(T) 



This estimate can be improved in two ways 
(i) Gamow Peak Asymmetry 



\ ' ( 




4V3 \ 




s define 





(3.22) 



I = Soe-'y°^—yoF{yo) , (3.23) 

where F(yo) can be evaluated by a term by term integration of the non-gaussian 
terms in the series expression of 0; for t/o >> which is typically the case for 
stellar plasmas, this correction can be written as a polynomial series in 1/yo 

F{yo) = l + ^ + ^ + ...^l + -^ + -^ + ... . (3.24) 
yo Vo 36yo 2592|/^ 
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(ii) Energy dependence of S 

As I/O depends on T, one usually writes 5* as a (truncated) Taylor series of starting 
point E = Q 



S{yT) = Y,s^y-T- 

n=0 

S(")(0) 



(3.25) 



where the parameters s„ = — should be fitted from experimental data; it 
is worthwhile to stress that, in general, S'(O) should not be understood as the 
zero energy limit of S{E), but rather as the zero-th order term in the polynomial 
approximation of the S factor in the energy range covered by the data, since the 
low energy behavior could show very different trends. 

Using Equation ()3.24|) and ()3.25p one gets 



dyS{yT)e 



n=0 



where 



(3.26) 



(3.27) 



Truncating the series at m = 2 and picking up only the 1/yo correction for F{yo), one 
derives 



{av) 



where 

SeffiEoiT)) 



So 1 + 



5T 



-Seffe T 



3T 



36Eo 

Finally the rate can be written as 



+ 51^0 1 + 



35T 
36En 



+ 52^0 1 



89r 
36En 



{av) (T) = A exp{-BT-^/'') (l + E ^n^"^') > 

^ n=l ^ 



(3.28) 



(3.29) 



(3.30) 



which is the analytical approximation often used in the literature. 

However, a better, though semi-analytical, approximation of the rate with respect 



to the Equations 

IniT) 



can be obtained as follows. We rewrite the /„(T) as 



I All 



where 



a„(T) = 



e-''''y^\ 



-yo 



(3.31) 



(3.32) 



Note that ao(T) coincides with the Fi^yo) function introduced to describe the Gamow 
peak asymmetry. This expression can be numerically evaluated in the relevant range 
and then fitted assuming the ansatz 



yo 



in) 



yl 



(3.33) 



Nuclear Reaction Network for Primordial Nucleosynthesis 



24 



somewhat reminiscent of the series expansion in the parameter y^^. The assumption 
?/o >> 1 is not always justified for BBN, nevertheless the previous ansatz reveals to be 
a good choice, and with the exclusion of the higher temperatures, the fitting errors are 
always less than few percent, typically less than 0.1% for n = and n = 1, which are 
the most relevant terms. The final form of the rate is then 



(-) = 1 /^^e-^ t SnEma^m , (3.34) 

usually reliable if one cuts at m < 3. This form, that was adopted in some cases in 
our previous work [201, is a good compromise between accuracy and physical insight. In 
the present analysis, however, we numerically integrate the (fitted) S{E) factor, often a 
polynomial function, in all the rate evaluations. 

Notice that in some cases (e.g. ^He + ^He 7 + "^Be) S{E) is better described by 
adding an exponential decreasing term. If 

52(£^) =^2(0)exp[-aE] , (3.35) 

then 

^2(0) 



/(2) = -jr- dEexp 




(3.36) 



from which, defining = T/(l + aT), we get 



T Att 

J = 5(0)^e-3^<y-yo , (3.37) 

where yo = {yc/'^)^^'^ and yc = Ec/Ta- We finally get in this case for the rate 

{av)2{T)=A,T~-"\Taf'eM-B,T~''^] , (3.38) 
with A2 and B2 to be (usually) fitted using experimental data. 

3.2.3. Resonant Reactions. The cr(-E) for an isolated (single-level) resonance is 
described by a Breit-Wigner formula 

- 2,^,E (E - Eny + {T/2y ' 
where uj is the statistical factor: 

(2J+l)(l+^„.) 

(2J„ + l)(2Ji + l) ^ ' 

while Fi is the partial width (PW) of the Compound Nucleus (CN) formation from the 
ingoing channel, and r2 the CN decay PW into the outgoing channel. Finally F and Er 
are respectively the CN total width and resonance energy. 

For F << En {narrow resonance) the integral ()3.14|) can be rewritten as 

M^(^j (^7)/?e-^ , (3.41) 
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where {uj'^)r = uVaVij/V and u is the statistical factor defined in (j3.4(J|) . This expression 
is a useful guide in fitting rates with narrow resonances. Note that no simple functional 
forms exists for the broad resonance case even if formulas like 

Co 



{av) = exp 



TV3 



[Eo{T)-En] + iT/2)' 



EAT^^' , (3.42) 



are sometimes used. 

Usually, for both neutron and charged-particle induced reactions, in presence 
of narrow resonances we added a lorentzian peak [1 + {E - Er)^ /{O.bVf]-^ to the 
polynomial component, while in the case of a single, broad resonance we multiplied 
the polynomial term by a similar lorentzian shape. In both cases, the integration 
was performed numerically and the fitting form usually chosen by looking at the past 
literature. 

Other cases are also possible (subthreshold or overlapping resonances), but they 
are of minor interest for our purposes and for further details the interested reader is 
addressed to the speciahzed literature. 

3.2.4- Screening and Thermal Effects. When treating nuclear processes in 
astrophysical/cosmological environments, one should carefully take into account the 
different physical conditions existing there with respect to that explored in a typical 
laboratory experiment. For example, the nuclear reactions are experimentally studied by 
accelerating partially or completely ionized species a on atomic or molecular (electrically 
neutral) targets h. As a consequence, the projectile feels a screened Coulomb potential 
because of the presence of the electronic clouds, that can be parameterized e.g. by a 
Yukawa-type potential 

0™(r) = ^e'^/^- , (3.43) 
r 

where Rah is an atomic/molecular effective radius depending on the projectile and 
(mainly) on the target properties. 

The standard phenomenological approach to this problem consists in introducing 
the screening energy parameter Ue = ZaZ^e^ / Rat, that represents an effective lowering 
of the Coulomb barrier seen by the projectiles, such that 

exp — . (3.44) 



As Ue are of the typical order of magnitude of the atomic energies, these corrections 
should be applied to correct the measured cross section for ii^ ~ 10 keV or lower, and 
thus only marginally affect the physics of BBN. Nevertheless there are cases, as for the 
+ p + or the ''Li + p ^He + ^He reactions, where this correction has to 

be properly taken into account"^. 

We note in fact that the value deduced for some reactions in deuterated metals seems to be one or 
even two orders of magnitude greater (up to 700 — 800 eV) than naive upper limits. See e.g. for 
an analysis of this phenomenon with several metals for the ^H(d,p)'^H reaction. 
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During BBN all matter is fully ionized and nuclear reactions take place in 
a surrounding plasma. This means that all cross sections deduced by laboratory 
experiments should be corrected for the effect of target-nuclei thermalization and 
screening. In fact, when temperatures are comparable to (or higher than) the excitation 
energies of the nuclei, all the accessible states are soon populated in a statistical way. 
Since the experimental determinations refer instead to the ground states, one has to 
calculate the needed corrections via theoretical models. For our purposes, we adopt the 
model used in ^2], based on the assumptions of local thermodynamic equilibrium of the 
plasma, and the use of the Hauser-Feshbach model to calculate the ratio between the 
reaction rate for the target in the n-th excited state and the one for the ground state 
nucleus. A further change should be also applied to the forward/reverse rate relation 
(see [Appendix E| ) where the statistical factor should be changed as follows 

^?,^(2J° + l)^5:(2Jr + l)exp(-|) , (3.45) 

if e" and are the energy and angular momentum of the n-th excited state of nuclide 
i. Useful fitting formulas can be found in Actually, the light nuclides have few 

(or not at all) excited states, so that for such corrections to be of some importance, 
(say ~1%) not only high temperatures, but also high mass numbers A are needed. This 
explains why, for the BBN purposes, the inclusion of this effect only marginally (i.e. 
< 1%) changes the standard predictions of the Lithium isotopes. 

We close this Section with some considerations on the screening effects for the 
nuclear reactions in the BBN plasma. The treatment of this problem is simplified since 
the plasma is non-degenerate and weakly coupled, i.e. the Debye screening length Rjj 
is much greater than the mean separation between particles. When the thermal 
pairs dominate the charged component of the plasma the enhancement factor / of the 
thermonuclear reaction rates has been calculated in [02] 

/ = exp (^^) = exp(1.28 x lO'^ZiZs), ^ « 1 , (3.46) 

so that the effect of screening is of the order of 0.1% at the highest temperatures involved 
in the BBN. Even smaller corrections are expected at lower temperatures (T < 0.1 
MeV), when the nuclear chain becomes effective and almost all pairs have already 
annihilated. In this case / can be calculated using the standard Debye-Hiickel theory, 
which gives 



with ( = + Zj)Xj ~ 2Xp + 6X4He, since during the late BBN phase protons and 

^He nuclei are largely dominant. Therefore 

/ = exp ^ exp(l.lxl0-V^10i0CZiZ2), ^ » 1 , (3.48) 

that is lower than what is obtained from ()3.46|) . as expected because of the very low 
baryon density. 
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3. 3. Error Estimates for Nuclear rates 

Unless it is possible to rely on a specific theoretical model, the S function for 
various processes*, both in magnitude and energy behavior have to be deduced from 
experimental data. This would not be a problem at all if were not for the fact that often, 
for a single reaction, it is necessary to combine several data sets affected by different 
statistical and normalization errors, and which usually cover different energy ranges with 
only a partial overlap. In our discussion we will use the following notation. With Ei^, 
Si^, aii^, and we denote, respectively the (center of mass) energy, S factor, statistical 
and (relative) normalization uncertainties of the i-th data point of the k-th data set. 
Whenever only a total error 0"*°* determination is available for a certain experiment, 
that error is used instead of ajj,, and an upper limit on the scale error is estimated as 
max[<7^.J. 

In many cases the different data sets do not agree within the quoted uncertainties, 
perhaps for the presence of scale/total errors larger than quoted. There are two possible 
solution to this problem. In one case there may be plausible arguments suggesting 
that one or more than one data set are seriously affected by systematics, which are 
however difficult to quantify. In this case a conservative approach simply consists in 
removing that particular data from the analysis, provided this would not severely limit 
our knowledge on the corresponding reaction rate, i.e. if there is no other information 
on S in the same relevant energy range apart from those data. On the other hand, 
lacking any hint to decide whether a particular experimental result gives biased results 
or not, there are basically no other choices than keeping all experimental information, 
but conservatively assume that data for all experiment may be still affected by an overall, 
albeit unknown, normalization factor. 

This approach simply means that, differently than in a recent analysis performed 
in ^H], we do choose not to fix the mean normalization for the data sets to unity, but 
rather introduce a normalization parameter Uk, which eventually have to be included 
among the fitting parameters, while we still keep the quoted as the estimate of the 
corresponding relative uncertainties. 

This procedure may seems quite arbitrary. Indeed, if we had no information at 
all on the possible energy behavior of a definite S factor, we may potentially introduce 
a large bias in data analysis. We do claim however that this is not the case at hand, 
since at least qualitatively we do have such a prior knowledge, namely that the S factor, 
apart from resonance contribution, is a smooth function of energy, for which a low order 
polynomial fit is a rather good approximationjj. We are therefore confident that our 
method can give quite a reasonable fit of S, as well as of the corresponding uncertainty. 

Other approaches have been used in the literature. For example, in j^l] a Monte 

* In the following we consider the S function, but all results also hold for the R function used for 
neutron induced reactions. 

jj We think that a theoretical prejudice is in any case entering in all possible data analysis of S factor 
and, we may say, of any physical observable. For an intriguing discussion on this controversial issue see 
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Carlo method for direct incorporation of nuclear data in the BBN code is adopted, 
which has the benefit of simply skipping the intermediate fitting step and provide a 
straightforward method for the inclusion of new data. Unluckily, this method also 
presents several drawbacks. For reactions with very little experimental information, 
the method is of very limited application. Moreover it may provide unphysical results 
in data-poor energy regions, and also underestimate the errors in intervals where no 
overlap among data sets exists. See also [12] for a detailed discussion on these and 
related issues. 

Very recently another detailed treatment of data for the BBN main reactions has 
been considered in which makes use of the R-matrix fitting technique. The authors 
stress the better parameterization of the cross sections obtained through this analysis. 
While this may be true in several cases, it is difficult in our opinion to claim any 
statement of generality, since all models rely in fact on some assumptions. For example, 
for non-resonant reactions the R-matrix model mimics the smooth behavior with some 
high-energy resonance tail, that actually is not physical, since a direct interaction 
contribution well described in the Born approximation would be more appropriate. In 
this case, the R-matrix approach could require more parameters than a (low-order) 
polynomial fit. Moreover, the R-matrix is known to work better in the extrapolations to 
very low energies often needed in stellar astrophysics: however this is of poor interest 
for BBN, where the reliabilty of the interpolation between the data over a broad energy 
range is crucial, and that can be equally (if not better) provided by simpler fits. 

Their error analysis, aimed to give a statistical meaning to the fitted quantities, 
makes a somehow arbitrary distinction between the case of underestimated statistical 
errors, which are then corrected with the standard protocol of infiating the x^, and that 
of discrepant normalization, where instead another procedure is used. 

Another issue which should be addressed is how to properly account for the existing 
correlation among data obtained by the same experiment (see e.g. jHl] and References 
therein). Again, this point is not as easy as one may think. One possible approach, as 
the one used in jl^, is to directly include the estimated normalization uncertainties in 
the covariance matrix which is used to construct the function. However, as discussed 
at length in jHH] , this does not seem a fully satisfactory choice, because of some subtleties 
related to standard error propagation. In fact, this procedure produce a negative bias, 
with the usual definition of bias in statistics, which systematically underestimates the 
physical quantities, in particular in presence of relatively large normalization errors and 
large numbers of data. Moreover it typically gives more weight to data sets affected by 
larger systematics, provided they have instead a smaller statistical error. We address 
the reader to jHSI for an extensive discussion on these issues. 

Our fitting method can be summarized as follows. We consider the function 




(3.49) 
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where 

Xstat — , ,2^2 ' yO.OU) 

Xnorm ^ ^ 2 ' (3.51) 

The a;fc are as mentioned free parameters of the fit, giving the renormahzation constant 
for the k-th. data set, while Sth{E) is a function of energy, typically a polynomial, which 
is chosen as the 5" fitting expression. It depends on some parameters a; to be determined 
by standard minimization procedures, along with the Uk- For the choice of the fitting 
function, usually done by looking at the pre-existing (theoretical and experimental) 
literature, see also section 13.2.31 and [Appendix D The sensitivity to the fitting form 
was checked by varying the number of free parameters (e.g., the polynomial degree) until 
one finds a "plateau" of the minima of the reduced the model space around some 
configuration of parameters. The model lying in this region with the lowest number of ai 
was finally chosen, in order to minimize unwanted fluctuations between the data. Notice 
that correlations among the data of the same experiment are taken into account in two 
ways. First of all, all data of some data set share the same value of Uk- Secondly, we 
add a penalty factor to the x^, given by Xnorm^ that disfavors a choice for Uk — l greater 
than the estimated normalization or total error ek- The fitting parameters coming from 
the minimization procedure are expected not to suffer of the bias problem discussed 
above, as it will be clear from its explicit application to several reaction rates in the 
next Sections. Note that our method is the most natural generalization of the unbiased 
one presented in [0^] (see Equation (3)), and the one actually suggested to be used (|SS1, 
Section 4). Notice that there are as many terms in Xnorm norm as there are in xltat- This 
takes into account that both the normalization and the shape uncertainty (in principle, 
on an equal footing) do contribute to the overall indetermination. 

Let us now discuss our estimate of rate uncertainties. From results of the 
previous Section, we recall that, through S{E) or R{E), the rate is a function of the 
parameters ai whose best estimate di and covariance matrix cov{ai,aj) are obtained 
via x^ minimization. The rate / is then obtained by numerical integration of S{E) 
(or R{E)) convoluted with the appropriate Boltzmann/Gamow Kernel (see Equations 

and srn^ ) 

POO 

f{T)= dEK{E,T)S{E,a) , (3.52) 
Jo 

while its (squared) error 6f trough the standard error propagation as 

if = r ,E' K(E\T) r dE K(E, T) E 

JO Jo Ti oaj 



cov{ai,aj) , (3.53) 

,j oai a oaj i 

thus fully including the correlations among the fitted parameters. As we said, a 
dominant role could be played by the systematic error. The systematic discrepancy 
between the data is shown by values of the reduced chi squared significantly greater 
than 1. Note that its value has been also partially increased because of our choice to 
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add x\ormi while on the other hand further parameters uj^ have to be estimated from 
the data thus reducing the number of degrees of freedom. Since in our approach the 
bulk of the systematic uncertainty related to the different normalization of the data has 
been taken into account, any residual discrepancy can be considered as due to some 
unidentified/underestimated source of error in one or several experiments. In this case 
we simply use the standard prescription of inflating the estimated error by the factor 
^J~)^ (see for example [3]). For illustrative purposes, for each reaction we will also quote 
the quantity e defined as 

,,^Ef^^K^ . (3.54) 

the sum being on the K different data sets, and the weights Wk chosen as Wk = 
{Xk/Nk)~^ , where xl is the contribution of the n-th data set (with A^^, data) to the 
X^. With this choice we assign a greater weight to the renormalization factors Uk of 
the data sets closer to the fitted function, which indeed is characterized by a lower 
value of Xk/^k- The value of e represents the typical renormalization of data, and is a 
qualitative estimator of the scale disagreement among several data set, but it should not 
be confused with a discrepancy error (that is already taken into account in our approach 
and already discussed) or a scale error, that we are now going to define. 

The overall scale error used in the analysis, was chosen to be equal to the lowest 
experimentally determined for that reaction. As it is by definition independent on 
E, and given the equation ()H.52|1 . it applies both to the S-factor and to the rate. It 
was added in quadrature to the statistical error in the fits presented in [Appendix E 



and in our following numerical analysis. The rationale under our choice for it is that 
combining several data set one should always expect a better estimate of the 5* factor, 
with a smaller error. Indeed this is not the case for data sets showing some serious 
discrepancy but, as we said already, this source of systematic error is separately taken 
into account. A different approach has been adopted in [TD . Apart from suffering in 
some cases of a negative bias on the best value determination, as we already emphasized, 
which typically leads to overestimate the discrepancy error as defined in Equation (22) 
of , in this analysis an intrinsic normalization error is introduced (see Equation (23) 
again in [12]). The latter represents an average of the normalization error, which is 
always greater than the best determination of the scale error. This typically results in 
larger estimates of the errors. 

A comparison with previous compilations generally shows a decrease in the error 
estimates (when available), and a shift in the best value of the rates for some reactions, 
typically of few percent. For the pre-NACRE era, this is likely due to significant 
differences in the datasets included, or to the use of analytical method instead of 
numerical one in the integration of the rates, and will not be discussed further. On 
the contrary, when a comparison is possible, our results essentially agree with the all- 
purpose NACRE compilation, though the inclusion of new data, the focusing on the 
BBN energy regime and the more "statistically motivated" method we used instead of 
the "minimum-maximum" rate given in jTH] naturally explain the lower uncertainties 
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found in our work. Qualitatively similar conclusions were found in (21] , and the residual 
quantitative differences with respect to our compilation can be fairly attributed to the 
slightly different databases and the different regression protocol used. 



3.4- Leading processes 

In this Section and in the following one we discuss in details all leading reactions in the 
BBN network, as well as a selection of those reactions which, though presently playing 
a sub-leading role, are affected by uncertainties large enough so that they still may 
contribute to some extent to the eventual nuclide abundances, or that present some 
historical or peculiar importance that make them worth to be discussed. The leading 
or sub-leading role of a reaction was firstly established graphically, by looking at the 
temperature behavior of their contribution to the right hand side of the corresponding 
Boltzmann equation for Xj (see [Appendix B and [14j for a similar approach), and then 



checked numerically. In particular for these two classes of reactions we report our results 
on adopted best estimates and uncertainties. A complete list of the full BBN network 
used in our study is presented in [Appendix A| For the sake of brevity, we do not report 
in this work all the statistical quantities entering each of the reaction analysis (see e.g. 
equations ()3.50|1 . ()3.51|) . ()3.54|1 ). The interested reader can obviously require further 
information by addressing a mail to one of the authors. 



3.4-1 ■ Reaction pn7.- p + n ^ ^ + . 

It is the main channel for production and, indirectly, of all other nuclides. There 
is a lack of experimental information about this non-resonant process, at least in the 
energetic range relevant for BBN. Except for the measurements of the thermal neutron 
capture cross section ath (see in particular [HE]), the only low-energy determinations 
are the ones in |E7j and jnH|- One can also rely on some deuton photo-disintegration 
data at energy near the threshold, as the historical ones of [Hll HOj (see also the review 
[7T]). that can be related to the direct process through the detailed balance principle. 
Recently, new photo-destruction measurements have been performed , that allow for 
an evaluation of this crucial reaction rate with an uncertainty smaller than 6% in the 
relevant range for the BBN, thus comparable with the estimate of but only based 
on experimental data. 

It is worth pointing out that this process has been theoretically studied well enough 
to allow for even better results. In fact, whenever a comparison has been possible, 
experiments have shown a substantial agreement with theoretical calculations j73j . 

An estimate of the rate based on the calculations of [Tl] and [73], normalized to 
ath = 0.332±0.002 6 (as measured by jZO]) was already presented in pUj . 
In J3] and in the previous BBN analysis the calculation of [77j has been used, which 
took into account the thermal capture evaluation of [ZHj {o'th = 0.3326±0.00076), on 
the high energy data of [79J and on the ones of deuton photo-disintegration, with a 7% 
total error on the rate obtained by adding in quadrature the evaluation errors (< 5%), 
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the fitting ones (5%) and the numerical integration approximations (2%). Note that in 
[Ti] it is erroneously quoted 0.006 6 instead of 0.00076 for the uncertainty on this cross 
section evaluation. 

In the present work, we evaluate the rate on the basis of the (few) available direct 
and inverse low-energy data and, in particular, the theoretical calculation of Rupak jHOl- 
As in the previous work of Chen and Savage [HI] , a pionless effective field theory is used, 
but the calculation is pushed to the next order, thus lowering the relative uncertainty 
to less than 1%. In this approach, since the relevant energy scale is much lower than 
the pion mass (~ 140 MeV), it is meaningful to describe the strong interaction among 
nucleons without explicitly introducing the pion degrees of freedom, and using effective 
four-nucleons local operators, while the electromagnetic coupling is obtained via the 
gauge principle. It can be shown that, at the energies relevant for BBN, the transition 
amplitude for the pn-^ process is dominated by the Mly (iso-vectorial magnetic dipole 
transition) and Ely contributions (iso-vectorial electric dipole transition), respectively 
at lower and higher energies; the Mly amplitude was calculated up to the next-to-next- 
to-leading order {N'^LO) and normalized to the thermal neutron capture cross section 
ath taken from |82]; the Ely amplitude was instead computed up to the N^LO order 
and normalized to the nucleon-nucleon scattering data and using the deuton photo- 
destruction measurements (in the 2.6 -^ 7.3 MeV energy range of the 7). When sub- 
leading effects are neglected, one obtains the effective range theory standard results. 




10 10 ^ 10 ' 1 

E (MeV) 



Figure 10. The low energy data and the fit of the auxiliary function R{E) for the p 
+ n ^ 7 + reaction. The solid line is the overall best fit, the dashed curve is based 
only on experimental data. The theoretical points taken from [SO] have a negligible 
error bar on this scale. The error bars for the experimental points also include the 
propagated uncertainty on the energies. 
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The fit of tlie R{E) function is almost completely determined by the Rupak results. 
The reaction rate was calculated analytically according to Equation ()3.18p . Linear error 
propagation leads to an estimate of uncertainty < 0.4% for Tg < 1.5. 

Since the database was limited io E <1 MeV, we expect that our rate will disagree 
with the previous estimates at high temperature; actually, starting from Tg ~ 2, with 
Tg = T/lO^/r, the difference with respect to the rate published in ^| has almost a 
linear growth; the inclusion of higher energy data in an auxiliary fit (~ 3 ^ 6 MeV) 
allows to check that 0.1% variations appear from Tg ~ 2. For this reason we choose to 
use our rate in the range Tg < 1.5, where it is certainly more accurate, and matching it 
with the rate evaluated in f^j, that is still used at higher temperatures Tg > 1.5. The 
overall uncertainty has been conservatively estimated as 1.2% for Tg < 1.5, and includes 
the theoretical (< 1%), the fitting (< 0.18%), statistical (< 0.4%) and normalization 
errors (~ 0.48%). The latter is due to the disagreement in the thermal capture cross 
sections |66j. The uncertainty grows to ~ 8% for Tg > 1.5, being based on the error 
budget of Pl] and the matching and normalization errors. For comparison, the analysis 
performed in jHHl gives a sample variance error of 4.45%, while the uncertainty quoted 
in [inj is 2.5%. All the available rates agree at about 1 a or better within the quoted 
errors. 

3.4.2. Reaction dp7.- "^H + p ^ 'j + ^He . 

It is one of the main responsible of the ^He synthesis both at high and low 
temperatures. When the deuterium formation channel is opening, this role is played 
instead by the strong reaction ddn., + ^ H ^ n + ^He (see Figure lB3|) . The dp^ is 
a direct capture reaction, for which it exists quite a recent theoretical model [HI]. The 
data sets we consider are [TH], [HS]-|n3]- Some discrepancies exist for the lower energies 
data, but in [90j the authors stress on the presence of a systematic error in their previous 
data (see e.g. [En]), with an upward correction that reduces the compatibility problem 
with the older data reported in [22]. Moreover, the latter data set by Griffiths is likely 
to be affected by a 10%-15% normalization error due to the wrong stopping powers used 
for the heavy-ice targets IHlinSl, so that the disagreement is now reduced with respect 
to the first claims. Finally, the recent data from the LUNA collaboration ^Hj allow to 
firmly establish the low energy behavior of the astrophysical factor: for example, the 
older analysis performed in |SS1, not including them and with a less accurate treatment 
of the other experiments systematics gave a sample variance error of 13.2%, while our 
present study suggests an overall uncertainty for the rate less than 3%. In the recent 
compilation ^H] the estimated error is close to 7%, but this number is dominated by the 
normalization error and so it is likely to be overestimated (see our discussion in Section 
3.3). We get a data normalization spread of e ~ 1%. As in our previous work |2IIj, the 
fit of all the data available has been performed with a cubic polynomial, and both the 
rate and uncertainty were obtained by numerical integrations, in fair agreement with 
the semi-analytical formula introduced in Section 3.2.2. 
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Figure 11. The data and fit for tlie astrophysical factor of + p ^ 7 + ^He. Tlie 
fit of S{E) is with a cubic polynomial. 

3.4.3. Reaction ddn; '^H + '^H^n + ^He . 

At Tg ~ 1, almost all '^He is produced through this typical Direct Interaction 
channel (see Figure lB3|) . The energy range of main interest for the BBN is 0.01 2.5 
MeV, with a particular relevance of the values E < 0.4 MeV. Apart from some windows, 
the experimental situation for both the ddn and ddp (see later) reactions is quite good, 
if compared with other processes: they are strong interactions, their Coulomb barrier is 
low, and their relevance as thermonuclear fusion energy routes made their study quite 
appealing. Nevertheless their importance for BBN requires further efforts. We used the 
data |M]- |100j : note that the data of |101j . used in ^H], are only a preliminary report of 
the measurements quoted in [^Hl, so they are not an independent data set and should 
be excludedtf. The data of [HS] quite accurately fix the high energy behavior of this 
astrophysical factor, and the data in [OHl allow to quote a normalization of the S factor 
(and then of the rate) of 1.3%. As many data are available, and they are generally 
affected by a small normalization or overall uncertainties (at the level of few percent), 
even relatively low discrepancies on the scale and/or the behavior of the S{E) will show 
up in the value of the reduced x^, which is of the order of 10 or greater at its minimum. 
Our method gives an overall uncertainty at the level of 1.3%, while for comparison the 
analysis performed in jHS] gives a sample variance error of 3.1%. It should be pointed 
out that, while the uncertainty expected on the basis of the best determinations should 
be even lower, when combining several data sets the experimental situation suggests the 

ft Incidentally we note that there is a misprint in 19 , as the data of jl()2j referring to the tdn reaction 
are erroneously quoted. 
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Figure 12. The data and fit for the S factor of the + ^ n + ^He reaction. 



presence of systematics. Indeed some experiments are likely to have underestimated the 
errors or could be affected by some unknown bias. In this situation, every regression 
method makes poor sense, and we continue to follow our approach just for consistency. 
As a confirmation of how difficult the analysis is in this case, we observe that the 
error budget estimated in a completely different approach in ^Hl is greater for the ddp 
6.9%) than for the ddni^^ 5.5%) reaction, despite of the fact that the quality and 
the source of data are comparable, when not better and even more abundant for the 
first reaction. 

As this reaction and the conjugate ddp strongly influence the error budget, 
so that it affects many of the nuclides predictions, we strongly recommend a new 
experimental campaign to determine accurately (say at the 1% level) both the magnitude 
and shape of the S{E) in the wide range useful for BBN studies (E< 2.5 MeV), and in 
particular for 0.2 < E <1 MeV, where actually no data presently exist. 

344. Reaction ddp.- + ^ p + . 

This process is the leading source of direct tritium synthesis (see Figure lB2|l . The 
discussion of its experimental situation closely follows our analysis of the previous 
reaction ddn, apart from the fact that very low energy data are available |lUUj . which 
have to be corrected for screening effects in order not to bias the result. The data we 
used for this reaction are contained in the same References quoted for the ddn, as these 
processes are usually studied together. Since these two reactions are non- resonant, in 
both polynomial fit for S{E) is used, while the rates and errors were obtained by 

numerical integration. The uncertainty in the relevant temperature range is less than 
1%. 
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Figure 13. The data and fit for the S factor of the + ^ p + ■^H reaction. 

3.4.5. Reaction tdn; + ^ n + "^He . 

The tdn is the fundamental channel for ^He synthesis during BBN. Many data 
are available for this process, also because it is the most promising process for the 
thermonuclear fusion research (low Coulomb barrier and high Q- value). Actually it 
is the energy production mechanisms considered for the next generation controlled 
fusion reactor, ITER jin3j . A broad resonance in ^He, having Eji = 0.064 MeV and 
r = 0.076±0.012 MeV is superimposed to the non-resonant background. This allows 
to use a Breit-Wigner shape to extrapolate with some accuracy the 5* factor also to 
energies below the range covered by the data. The data sets we used are |l()4j - PTT] . 
with the Conner data jl05j excluded at energies greater than 240 keV, where the isotropy 
assumption for cr[E) fails. The rate and error estimate were both calculated numerically. 
Our error estimate is ~ 0.6%, and e ~ 0.02 (for a. xl — 1-4), while for comparison the 
analysis performed in [Sn] gives a sample variance error of ~ 4%, while the discrepancy 
error estimated in [THj is about of 2.2%. 

3.4.6. Reaction heSdp.- ^He ^ p + ^He . 

Like the conjugate reaction tdn, it is dominated by a broad ^He resonance, in this 
case at ~ 0.2 MeV. There is quite a satisfactory experimental knowledge of this process, 
even if the dispersion of the data grows for E < 0.2 MeV, introducing some uncertainty 
on the resonance parameters. It is the second route to ^He production after the tdn 
reaction, but its indetermination mainly plays a role in the ^He and ^Li yields, as it 
essentially controls the burning of ^He. The data considered in our analysis are the 
ones reported in |112j - P^K] : the two recent data sets |124| 112,^^ . firstly included in this 
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Figure 15. The S{E) data and fit for the ^He + ^ p + ''He reaction 



compilation, once corrected for the screening effect allow for a better coverage of the 
low-energy region. Our new rate agrees within the errors with the rate published in 
|14j . where the total error is estimated to be 8% (1 a). For comparison, the analysis 
performed in [521 gives a sample variance error of 9.15%, while in a 7.3% result is 
quoted. We estimate a rate uncertainty of the order of 0.6%, while e ~ 0.03. 
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Figure 16. The S factor data and fit of the ^He + n ^ p 



3.4-7. Reaction heSnp; ^He + n ^ p + . 

This process keeps '^He and in equilibrium at (relatively) high temperatures. 
During BBN it plays an important role in determining the final abundances of ^He and 
''Li. In our regression we used the data reported in |126j - [T!?n] . by limiting the analysis 
to E < IMeV as in this allows to cover the most interesting range for the BBN by 
significantly reducing the numbers of parameters needed to have a good fit. Thanks to 
the new measurements of the reverse reaction cross section near the threshold reported 
in |13nj and to the accurate knowledge of the thermal capture rate |128j . both the 
normalization and the shape of the R{E) factor are well known, so that its contribution 
to the final error budget is low enough now not to require further measurements. In 
particular we find a statistical uncertainty of at most 0.2%, and a normalization spread 
of e ~ 0.002. 



3.4-8. Reaction be7np.- ^ Be + n ^ p + Li . 

It is the main responsible for ^Be/''Li conversion during the relevant phases of the 
BBN (in the final stage, the ^Be electron capture is actually more important). Only for 
very low energies, from 25 meV to 13.5 keV, data are available on the direct process 
|131j . while for higher energies one has to rely on the data for the reverse reaction 
|1321 1133] . To avoid the introduction of significant errors due to the use of the detailed 
balance conversion near the threshold, we considered indirect data with energy > 40 
keV only. As for the previous he3np reaction, we also restricted the analysis at energies 
E < 2 MeV, that are essentially the only ones relevant for BBN. However, differently 
than the non-resonant heSnp reaction, whose rate was obtained analytically, this process 
has a resonance at £^ ~ 0.32 MeV that suggests to use a fully numerical approach, both 
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Figure 17. The S factor data and fit of the ^Be + n ^ p + ^Li reaction. 

for the rate and the uncertainties. Note that despite of the few data set available, the 
Koeler data |131j fix the overall scale error to the ~ 2.1% level, thus making this process 
quite accurately known for the purpose of BBN studies. The statistical error is of the 
order of 0.7%. 

3.4.9. Reaction Q;he37.- ^He + ^He ^ 7 + '^Be . 

It is the dominant channel for direct ^Be production and, in the (relatively) high-r; 
Universe suggested by WMAP data, practically all ^Li synthesis is controlled by this 
reaction. Its importance for the prediction of the solar neutrino spectra has motivated 
several theoretical and experimental efforts in the past years to obtain a better estimate 
of the cross section. A constant plus a decreasing exponential times a polynomial 
was used to fit its non- resonant 5* factor, and the data used are |134j - [1^T] . where 
the data in |137j were renormalized by a factor 1.4 to correct for the Helium gas density. 
Our regression method estimates an error of less than 3%, mainly dominated by the 
scale uncertainty. Both the value of the reduced x^; xl = 2.1, and the high average 
renormalization we find from the fit, of the order of 20%, in agreement with the majority 
of the quoted errors, strongly suggest to undergo a new measurement campaign, to finally 
establish both the shape and the scale of this process at a few percent accuracy level. 
Of course the solar neutrino predictions would benefit of this new data, too. 

3.4.10. Reaction \i7pa: Li + p ^ ^He + '^He . 

This is the main process destroying ''Li during the BBN. The data used in our 
analysis are the ones quoted in jl42| - p"i6j . A self-consistent correction for the screening 
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Figure 18. The S factor fit and data for the ^He + '^He ^ 7 + ^Be reaction. 
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Figure 19. The data and fit for the S factor of the ^Li + p ^ '*He + ^He reaction. 

potential was also performed, whose effect is particularly significant for the low-energy 
data of Engstler et al |145l I146j . Our estimate for the overall error is of the order of 2% 
in the relevant temperature range, while e ~ 0.01. 



34.11. Reaction ad^i: ^He + ^Li . 

Even if it is a weak electric quadrupole transition, this reaction is important as 
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represents the main ^Li production reaction. It is experimentally studied only at E > 
1 MeV and around the 0.711 MeV resonance. The low energy range is only weakly 
constrained by an upper limit, while the theoretical estimates for the non-resonant rate 
compiled in show differences of 1 ^ 2 orders of magnitude and were used to establish 
upper and lower limits. These features suggest the introduction of a temperature- 
dependent asymmetric uncertainty. It would be useful to get new data at E < 0.6 MeV 
in order to establish a reliable estimate of the standard BBN production of ^Li; in fact 
one cannot presently rule out the possibility that a relevant fraction of the observed ^Li 
in PopII stars has primordial origin (see |147j and our discussion in Section 4.1.5). In 
view of the serious lack of data, we simply updated the code with the inclusion of the 
NACRE evaluation for rate and errors. 



34.12. Reaction h6phe3.- ^Li + p ^ ^He + ^He . 

It is the main ^Li destruction channel. Available data for the cross section are quite 
abundant and accurate 

We used the data |148j - [135] as well as those in References |143tll45j . already cited 
for the lilpa process, with the last one self-consistently corrected for the screening effect. 
The estimated error is lower than 2% and, quite negligibly contributing to the eventual 
accuracy of ^Li estimate, with respect to the one coming from the adry reaction. We 
estimate e ~ 0.09. 

3.5. Main sub-leading processes 

3.5.1. Reaction at'^: ^He + 'j + Li . 
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Figure 21. The data and fit for the S factor of the "He + 7 + Li reaction. 



This is the main channel for the direct synthesis of Li, and has been known since 
a long time as a crucial process for the ^Li predictions of the BBN (see fl]). However, 
especially thanks to the recent data of |154j . its rate is relatively well known and, within 
the uncertainties, it only contributes at the percent level to the ^Li yield. This is one 
of the main consequences of the present suggested range for cj^, that indeed points out 
the leading role of the aheS'y reaction as the main route to A=7 elements production. 
We parameterize the S factor as in the case of the aheS'-y reaction. The data used are 
|154j - jl57j . and our regression protocol gives xl — 0.71, with an overall uncertainty 
estimated at the level of 15% and e ~ 0.02. 



3.5.2. Reaction tp^y: '^H + p ^ ^ + '^He . 

This non-resonant process is the third one in order of importance involved in the 
^He synthesis, but its influence on the final error budget of all nuclides is marginal. The 
old rate for this reaction is that published in the compilation ^2]- Meanwhile, new data 
were taken |1581 1159j . and a new fitting of the astrophysical S factor is now available 

S{E) = So + SiE + S2E^ , (3.55) 

50 = 2.0±0. 2 keVmb , (3.56) 

51 = (1.6±0.4)xl0"^m6 , (3.57) 
^2 = (I.l±0.3)xl0-^m6/fce\/ . (3.58) 

We used these parameters to estimate the best rate, while its uncertainty is assumed to 
be A5'(0)/5'(0) ~ 10%. Notice that correlations between the parameters haven't been 
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Figure 22. The data and fit for the S factor for the ^Li + p ^ 7 + *Be reaction. 

published, so our error propagation method on this rate cannot be used without a full 
re-analysis of the data. 

3.5.3. Reaction h7p7.- Li + +^He + *He . 

This reaction mainly proceeds trough the resonant term at -EsBe = 19.86 MeV, 
but also has a non-negligible non-resonant contribution. Both were quite recently re- 
determined in The relative importance of this reaction as an additional burning 
channel for the ^Li was pointed in ^3] , where its contribution was estimated to change 
the destruction rate up to ~ 8% at Tg = 1, but this seems to be neglected in all recent 
studies. Notice that the new data collected in |16Uj move upwards the estimate of the 
non-resonant contribution by more than a factor 10, thus further increasing its role. 
Tough sub-leading, it is worthwhile to note that within the actual uncertainties and 
assuming the WMAP range for tUb, this process acquires a role comparable to that of 
^Li +p — > ^He + ^He and ^He+^H— > 7+^Li reactions in determining the final prediction 
of ^Li. The data sets used in our analysis are |16Uj - P^^ . We estimate an overall error 
of 3% and a large value of the normalization spread parameter, e ~ 0.14. 

3.5.4. Reaction he7na: ^ Be + n '^He + '^He . 

To our knowledge, evaluations for the rate of this reaction have only been published 
in P and [2Sj; without information on the sources of the data and error estimate. We 
did not find further analysis in the following compilations by Fowler et al . ^0]- The 
two data sets of the reverse process published in |163l I164j refer to center of mass 
energies of the direct one greater than 0.6 MeV, thus leaving a great uncertainty in the 
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BBN window. They seem to be roughly consistent with the old estimate of the rate, 
and a new one in view of so scarce data would make little sense. For this reason we 
adopted Wagoner's rate, assuming a factor 10 uncertainty, as he suggested as a typical 
conservative value. Within this allowed range, this reaction could play a non-negligible 
role in direct ^Be destruction, so it would be fruitful to have a new experimental 
determination. Apart from the role of unknown or little known '^Be resonances, it 
is however unlucky that the used extrapolation may underestimate the rate by more 
than one order of magnitude, as this process mainly proceeds through a p-wave. 

3.5.5. Reaction h7dna.- '^Li+'^H^n + "^He + ^He . 

The deuterium induced destruction processes of ^Li are dominated by this reaction, 
and a relatively recent analysis on this issue has been published in Reference |lfi5j . 
However, within a conservatively assigned 50% uncertainty, the new evaluation produces 
almost no change in the Lithium amount with respect to the compilation J2] 5 so it seems 
difficult that this process could play a role in reducing the ~ 3a discrepancy between 
predicted and observed ''Li yields. 

3.5.6. Reaction be7dpa.- '^Be + '^H ^ p + ^He + '^He . 

The rate for this reaction is taken from ^2], and is probably derived from the 
Parker estimate based on the data published in |l(i7j . It has been recently pointed 
out |168j that an enhancement of this process at low energies {Ecm < 0.5 MeV) of a 
factor 100 or larger could reduce the discrepancy between the predicted and observed 
values of ^Li. Even if this possibility cannot be ruled out on the basis of the available 
data, we stress that there are no compelling reasons, both theoretical and experimental, 
suggesting such an anomalous enhancement at low energies. New measurements will 
greatly help in clarifying this issue. 

4. Light nuclei abundances 

In this Section we first review the present data on primordial abundances of light 
nuclides. The aim is to provide a sufficiently self-contained description of the issue, 
mainly collecting the most recent results more than entering into the complexity of the 
measurement procedures and of all possible sources of systematics. For this reason it 
should not be considered as a complete discussion on such an involved topic. More 
details on the experimental techniques can be found in the original papers we do quote 
in the following. The second part of this Section is then devoted to compare the 
theoretical expectations, obtained by numerically solving the BBN system of equations, 
with experiments. Particular stress is given to quantify uncertainties, in particular those 
due to propagation of nuclear rate measurement errors. 

4.1. Experimental results 
4.1.1. . 
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The deuterium is known to be the nuchde that better allows for a precise 
determination of the value of rj. The most recent experimental results are discussed 
in |169j . where the authors deduce the estimate 

X.h/X, = (2.78t°:^^)-10-^ . (4.1) 

This value is the average of five measurements of DI/HI column ratio in different 
QSO absorption systems at high red-shift, obtained via deuterium isotopic shift from 
hydrogen. In particular this analysis takes advantage of the new detection towards 
Q1243+3047, based on an improved modeling of the continuum level, the Lya forest and 
the velocity structure of the absorption, as well as on a better relative flux calibration. 

In spite of the fact that the reduced is pretty large = 4.1, when using the set 
of QSO measurements considered in |l(i9j . so that unidentified systematics may be still 
affecting some of the data (correlations with the column density, local peculiar depletion, 
etc.), it is nevertheless interesting to note that the preferred range for X2H/Xp is in quite 
a good agreement with a standard BBN scenario with three neutrinos and WMAP result 
for ujb. We will come back to this point later. 

4.1.2. ^He . 

The ^He could provide other constraints on the BBN scenario, but here the 
experimental study is complicated by the fact that a relevant contribution is expected 
from stellar destruction and production, making difficult to extrapolate the primordial 
abundance of this nuclide. Apart from Solar System measurements, '^He is usually 
detected through the 8.665 GHz spin flip transition of '^He^, but data are limited 
to (chemically evolved) galactic HII regions and planetary nebulae (for this issue, 
see |17()j and References therein.) Assuming the simplest scenario of a typically 
limited stellar enrichment of ^He, the best observed system (S209 region) suggests 
Xa-fie/Xp < (1.1±0.2)-10~^, while a more conservative limit based on a larger sample 
gives Xsue/Xp < (1.9±0.6)-10"^. As discussed in |171j . at the moment one cannot 
really consider the '^He as a strong "baryometer" because of the uncertainty on the 
role of stellar processing, the limited range in metallicity explored, and also the weaker 
dependence of XsHe on rj with respect to X2h. On the other hand, if combined with 
the increasing precision of the BBN predictions, and using the value of rj as determined 
by present and forthcoming CMB data, these observations could clarify several issues 
related to the Galactic Chemical evolution. 

Finally, we note that the evidence collected till now for a sort of 'plateau' with 
respect to metallicity |17Uj and the order of magnitude of the observed abundance at 
least qualitatively support the idea of a common, primordial origin for the bulk of this 
nuclide abundance. 

4.1.3. ^He . 

The value of Yp is usually derived by extrapolating to zero metallicity measurements 
done in dwarf irregular and blue compact galaxies (BCGs), that are among the least 
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chemically evolved galaxies. The values obtained have a typically low statistical 
uncertainty (at the level of 0.002, i.e. less than 1%) because of the large number of data, 
but are affected by several systematical uncertainties, whose origin has been considered 
in many recent works (see e.g. References |172|ll73j ). The estimate given in |174j . based 
on a wide database and on a careful study of the chemical evolution of BCGs, is 

yp = 0.238±(0.002),t„t±(0.005),,, , (4.2) 

while in the recent work |175j the authors quote 

= 0.2421±(0.0021),t,t , (4.3) 

when the linear regression is performed vs. 0/H on a restricted, detailed-studied sample 
of galaxies. These observations are obtained with high signal-to-noise ratio, and a self- 
consistent approach on several He spectral lines was used in order to fit simultaneously 
as many parameters as possible. They also found that the typical size of the systematic 
errors is of the order of few percent, but the best estimate of the combined effect is at 
the 0.5% level. Slightly higher values (from 0.243 to 0.247) are found neglecting some 
corrections, using a wider data sample, performing regression versus N/H, or finally 
using different Equivalent Widths (EWs). 

Very recently, in |176j a new, detailed analysis of the budget of systematics in Yp was 
performed. With a different method and sample selection, this analysis of the same data 
of Izotov and Thuan gives Yp = 0.249±0.009, quite different indeed, even if consistent 
within 1 cr. If we adopt the conservative range quoted in |176j . (0.232 0.258), and 
viewing these bounds as a reasonable 2 a interval for Yp, we have (mimicking a "1 cr" 
error) 

Yp = 0.245±0.007 , (4.4) 

which indeed is the value we will use to compare our theoretical results. Even if a 
knowledge of Yp at the level of 3% may seem quite satisfactory, in order to use "^He 
to bound possible BBN scenarios, it would be quite important to improve the Yp 
measurement accuracy at the level of the quoted statistical uncertainty, which requires 
a clear understanding of main sources of systematics. 

We would briefly point out that other methods have been proposed to estimate 
Yp, that hopefully could help in understanding the role of unknown systematics in 
the previous results. For example, in \177\ and References therein, Yp is determined 
from indirect studies of Galactic Globular Clusters. The value they found Yp = 
0.250±(0.006),iai±(0.019) sys is still far from the required accuracy, but at the moment 
it is a promising contribution to clarify the debated issue of the primordial "^He 
determination. Finally in |178j a first attempt was made to determine Yp from CMB 
anisotropics data, but the 1 a confidence level 0.160 < Yp < 0.501 is still too broad. 

4.1.4. '^Li . 

The role of possible depletion mechanisms of the Spite plateau found in PopII 
dwarfs halo stars in the seminal paper |179j by Spite & Spite, is still under debate. 
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For example, the authors of |18(Jj found evidence for a positive dependence of Xtlj on 
metaUicity, estimating 

X7Li/Xp = 1.23±0.06l°;i , (4-5) 

at 95 % C.L., with uncertainty dominated by systematics effect. In other cases 
|18H 11821 I183j a negative trend and/or a dispersion of the measurements around 
the plateau greater than the experimental error was suggested. A huge number of 
mechanisms has been invoked to explain such a ^Li depletion, among which rotational 
mixing, diffusion, stellar winds, etc (see e.g. |184j and References therein for an 
introduction to the current status of PopII stellar models). However, according to 
|185j . these effects are second order, and a general consensus on the primordial origin 
of (at least) the bulk of the ^Li plateau has been achieved. The value quoted in |185j is 
the following 

X.Li/X, = 2.19^°:^^ . (4.6) 

An important step is the recent check that other objects like the Globular Cluster NGC 
6397 turn-off stars share essentially the same |185| 1186] , and even the troublesome 
observations in the globular cluster M92 reported in |187j seem to agree much better 
with other determinations in the reanalysis performed in jl88j with a new temperature 
calibration. 

At the moment, how to compare the extracted value of Xtl; with the one 
obtained from the standard BBN scenario is still quite a problematic issue. As the 
only traditionally observed line is the 6708 A doublet of the neutral Lithium (Li I) 
and the fraction of ionized Lithium (Li II) is estimated to be larger than 99.7%, 
a critical parameter for the extraction of A(Li)= log^gl'^Li/'^H) + 12 is the effective 
temperature of stellar atmospheres Teg. Other parameters as surface gravity, metallicity 
or microturbulence are much less important |186j . As Teg is only estimated indirectly, 
recently it was also pointed out the possible bias in extracting the plateau abundance 
from one spectral line only. Actually when the subordinate line at 6104 A is tentatively 
used, some disagreement is found, and values obtained for A (Li) could be high enough, 
A(Li)~ 2.5, to reduce the discrepancy with BBN predictions (further details can be 
found in |189j and References therein). 

As in previous studies, we do not use Xtli in our BBN likelihood analysis, waiting 
for a clearer understanding of possible systematics in stellar modeling. However, further 
efforts should be also devoted in trying to reduce the large theoretical uncertainty which 
is still affecting Xtli, mainly by lowering the uncertainty on the leading reactions, ^Be 
+ n ^ ^He + ^He, ^He + ^He ^ 7 + ^Be, and ^Be + ^H ^ p + ^He + "^He, which are 
responsible for a large fraction of the overall error on Xtlj. We will come back to this 
point in Section 4.2.3. 

4.1.5. ^Li . 

Among the remaining light nuclides, the most abundant (for the best range for rj) 
is ^Li, with Xau/Xp ~ 10"^"^. One has to conclude that, in the more viable scenarios. 
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there is no hope at the present for its detection. Actually, some positive measurements 
in halo stars at the level of ^Li/^Li^ 0.05 were obtained in the last decade |190j . but 
they reflect the post-primordial production of this nuclide in Cosmic Ray spallation 
nucleosynthesis. This does not mean that this issue is unrelated to BBN studies; on the 
contrary, as discussed e.g. in the study of the chemical evolution of the fragile 

isotopes of Li, Be and B could constraint the ''Li primordial yield, hopefully further 
clarifying the observational situation in the Spite Plateau. For example, the survival of 
such a (relatively) high fraction of ^Li in halo PopII stars severely limits the possible 
depletion factor for the more tightly bound ''Li. 

4.2. Theoretical predictions and uncertainties 

4.2.1. General considerations . 

The goal of a theoretical analysis of BBN is to obtain a reliable estimate of the 
nuclide abundances Xi and of their uncertainties, once the best values and errors on the 
standard physics inputs are known. These inputs are 

• Tn, the neutron lifetime; 

• Gn, the Newton gravitational constant; 

• f], the baryon to photon number density ratio; 

• the nuclear rates. 

The first two parameters are now known with a satisfactory accuracy ^ 

Tn = 885.7±0.8s , (4.7) 
Gn = 6.7087±0.0010- 10'^^ GeV'^ . (4.8) 

Until recently, BBN has been the only probe of rj, and the best fitting with the 
observed abundances was used to guess its value plus some insights on neutrino or non- 
standard physics. Thanks to the precise determinations of ujb obtained by the WMAP 
Collaboration results on CMB anisotropics, even rj is now independently known with 
an increasing precision, fixing the predictions of the standard scenario of BBN within 
a relatively narrow range. The value of the baryon parameter Ub = flbh'^ from CMB 
measurements is presently the following 

Ub = 0.023±0.001 . (4.9) 

As these first three sources of error are independent and uncorrelated, one can easily 
work out the propagated uncertainty on each Xi, by adding in quadrature the three 
contributions. 

The error matrix due to nuclear rate uncertainties are calculated as in our previous 
work [Hj 

4 = 7E[(^^(r^)-^^(r.)][(^.(r+)-x,(r^)] , (4.io) 
^ k 
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where the sum is over all reactions and with all other parameters set at their best values. 
The value is the average rate plus (minus) the estimated 1 a uncertainty. Thus 
the contribution to the uncertainty on Xj is given by 

a, = V^ , (4.11) 

and the correlation between Xi and Xj 



Being a generalization of the method described in |192j . this approach avoids the 
use of (time consuming) Monte Carlo simulations, while it does not require the existence 
of the functionals Xij{r]) = din Xi/dlnTk- This means that it can be applied even to the 
case of asymmetric or temperature dependent errors, which indeed is the case for several 
process rates. The reliability of the method is supported by the Central Limit Theorem 
(at least when many reactions contribute to the error almost on the same level) and the 
fact that there is typically a linear relation between the Xi and the parameters they 
depend on. Furthermore a direct comparison with Monte Carlo simulations results done 
in |192j shows that the method gives accurate results. We expect that this agreement 
is even improved for the typically smaller errors we found, and the allowance for the 
temperature-dependence of these quantities. 

Apart from the inputs we have discussed so far, a BBN scenario does depend on 
several parameters which set a specific theoretical framework. As a main example, the 
eventual abundances of (mainly) ^He is influenced by the fact that extra particle species 
contribute to the overall energy density at the MeV temperature scale. This increases 
the value of the Hubble parameter, so that neutron to proton ratio freezes out at larger 
temperatures, where the neutron fraction is higher. It is customary to parameterize 
their contribution in terms of the effective number of extra degrees of freedom, AN, 
defined as follows 

where p is the total energy density. Actually AN is a simple constant term for species 
which are fully relativistic at the BBN epoch, while it is a scale factor depending 
function otherwise: for example, for particles with mass much larger than MeV, AN 
increases proportionally to the scale factor. For the case of extra relativistic degrees 
of freedom decoupled from the electromagnetic plasma, the maximal allowed range 
0.232 < Yp < 0.258 quoted in [110] would constrain -1.14 < AN < +0.73 for 
Ub = 0.023. Since AYp ~ O.OISAX, a shift of 0.001 in would produce only a change 
of the order of 0.04 in the previous limits (see the Tables El and B] and [Appendix F ) . 

In addition to AN, light nuclei yields are in general sensitive to several parameters 
characterizing particle interactions or rather features of their distribution function. Most 
studied examples are the possible value of neutrino chemical potential |193j (for most 
recent bound on this parameter see for example \19^ 1^). or neutrino lifetime (for a 
review on this issue see [0]). 
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In what follows we will consider only the minimal standard scenario. Wc therefore 
assume no extra parameter at all, and so all theoretical predictions for ^H, "^He and ^Li 
abundances are determined by the value of rj. A careful study of this case seems to us 
of particular relevance since, as we mentioned already, the baryon to photon density is 
independently fixed by CMB anisotropy measurements, so there are no free parameters 
left. Comparing experimental result with theoretical expectations is therefore a clean 
test of consistency of the simplest BBN dynamics, based on our present understanding 
of fundamental interactions and of the cosmological model. It may also give useful hints 
to point out possible systematics in the experimental results reviewed in the previous 
Section. 

To compare the present sensitivity to of BBN with the one of WMAP, we will also 
study how this parameter is constrained by observed nuclide abundances, i.e. leaving 
its value as the only free parameter to be fitted, rather than using the WMAP value 
prior. To do this, or more generally to bound the values of any other set of parameters 
of the kind discussed so far, the standard procedure is to construct a likehhood function 

£{7j) a exp {-x\v)/^) (4.14) 

with 

X'iv) = UMV) - X:''^]W,,{r,)[Xj{r,) - ] , (4.15) 

where the proportionality constant can be obtained by requiring normalization to unity, 
and with Wij{ri) denoting the inverse covariancc matrix 

mAv) = [4- + <e.pS^J + 4-,o*.e.]-' , (4-16) 

where ai^exp is the uncertainty in experimental determination of nuclide abundance Xi 
and by crfj^other denote the propagated squared error matrix due to all other input 
parameter uncertainties (r^, Gn, etc.). 

4-2.2. Some remarks on rj . 

There are two different way to parameterize the baryon content of the Universe 
today, which can be used interchangeably, namely rj and cui,. It is interesting to review 
the simple relationship between these two quantities, in order to asses how uncertainties 
propagate from one parameter to the other. It is straightforward to get 
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This expression leads to 
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(4.18) 
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The uncertainties on the parameters (1^, Gat, ^70) generate a Arj/rj at the level of 
0.1% (errors added in quadrature), which in view of the present determination of ut, 
Aub/ub ~ 4%, is sub-leading and can be neglected in what follows. 

As a further remark we recall that the value of r] lowers from early times to the end 
of nucleosynthesis, because of the entropy transfer to photons during annihilation 
phase. Typically [22] the final value rjf is chosen as the input parameter, while the 
initial value, at T ~ 10 MeV, is then deduced as 

where we recall z = aT. A careful estimate of this relation should take into account the 
mass of the electron even at initial condition, so that entropy is not simply proportional 
to z^. Furthermore the value of photon temperature T get modified by the QED 
thermodynamical corrections to energy densities and pressures of and 7, as well 
as from the the small effects of the residual neutrino coupling to the electromagnetic 
plasma. Including all these one gets the estimate r]i ~ 2.73?7/, slightly different than the 
simplest result one would get neglecting these effects, and which is for example used in 
[22], Vi = 2.75r]f. Actually the corresponding corrections on the nuclide abundances are 
quite small. For Ub = 0.023 we estimate 

10^AX2h/Xp = 0.03, AYp = -0.6xl0-^ , IO^^AXtli/X^ = -0.07 , (4.21) 
In what follows f] always stands for rjf. 

4-2.3. Results with the WMAP prior on uob . 

As we mentioned, in the standard BBN scenario with no extra effective degrees of 
freedom or neutrino chemical potential, the light nuclide abundances are determined 
by the baryon density only. In Table [HI we report the values for the most relevant 
nucleosynthesis products, assuming the WMAP result, Ub = 0.023±0.001. We also 
show the corresponding uncertainties due to the propagated errors an on nuclear rates, 
as well as a^j^, the one induced by Wf, 1 cr error. As a general comment on these results we 
can say that our analysis of experimental rates, and in particular the adopted method 
for error estimates, led to a sensible reduction of an compared to our previous results in 
[20] • The reason for this is easily understood. In [20] we used a much more conservative 
approach, since we accounted for the fact that, for several rates, experimental results 
were often in disagreement by introducing a large maximal error defined similarly to the 
parameter e of Equation ()3.54p . 

As we already mentioned in the Introduction the value of Ub adopted in our analysis 
is the result given by the WMAP Collaboration combining their findings with those on 
CMB anisotropics of CBI and ACBAR experiments as well as with 2dFGRS data on 
power spectrum assuming a ACDM model. Slightly different values are obtained by 
including Lyman a forest data or introducing a running spectral index (see Tables 7 
and 8 in PJ). In particular the best value for Ub varies in the small range 0.022 0.024 
when combining different data sets, while the la uncertainty is always of the order of 
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nuclide i 


central value 






X2h/Xh (10-5) 


2.44 


±0.04 


+0.19 
-0.16 


XsHe/XH (10-5) 


1.01 


±0.03 


+0.02 
-0.03 




0.2481 


+0.0002 
-0.0001 


+0.0004 
-0.0004 


X.u/Xn (10-14) 


1.1 


±1.7 


±0.07 


Xru/Xn (10-10) 


4.9 


±0.4 


±0.4 



Table 3. The light nuchde abundances for ujb = 0.023 (second column). The (1 cr) 
uncertainties due to nuclear rates errors (an) and LUt (ffub) s-re also shown in the third 
and fourth columns, respectively. 



0.001, or slightly smaller. The impact of these different determinations of Ub on the nuclei 
abundance is shown in Table |3] where we report the best values for Uf, = 0.022, which 
is obtained combining WMAP±CBI±ACBAR in the running spectral index scenario 
PP, and for Uf, = 0.024, which is instead the preferred value of WMAP data only with 
a power law ACDM model. Notice that Deuterium and ^Li show appreciable changes, 
at the level of 10 %. We also note that the uncertainty due to nuclear rates, which in 
principle depends on the value of ui,, have very tiny variations and are still given by the 
values shown in the third column of Table 3. See also [Appendix F where we report a 



fit of both nuclei abundances and the covariance matrix for ^H, ^He and ^Li. 

Comparison with the analysis of ^H] shows a similar trend, though the uncertainties 
quoted in |^ are typically larger for the reasons discussed in Section As a main 
exception we mention the fact that we perfectly agree on the uncertainty quoted for Yp, 
since in this case (almost) all error is due to neutron lifetime statistical uncertainty. 

We now consider in more details our results for each nuclide. 

- 2H- 

The Deuterium abundance is in quite a good agreement with the experimental 
result quoted in Section 4.1.1. This is perhaps the main indication that indeed the 
standard picture of BBN is fully satisfactory. We know in fact that is extremely 
sensitive to the value of baryon density. Thus it is rather remarkable that the 
value found for Ub from a completely independent observable, such as the CMB 
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nuclide i 


ujb = 0.022 


Ub = 0.024 


X^h/Xh (10-5) 


2.63 


2.28 


X3He/XH (10-5) 


1.03 


0.98 




0.2477 


0.2485 


X.u/Xu (10-1^) 


1.1 


1.0 


Xru/Xu (10-10) 


4.5 


5.3 



Table 4. The light nuclide abundances for ujb = 0.022 (WMAP+ACBAR+CBI with 
a running spectral index) and ujt, = 0.024 (WMAP data only with a power law ACDM 
model) . 



rate 






49 


ddn 


37 


ddp 


14 



Table 5. The contribution of reaction rate errors to the total (squared) uncertainty 
on X2h/A"p 

anisotropies, leads to a yield very close (1 a) to measurements. 
In our findings the uncertainty due to nuclear rates is now reduced to 1.6%, mainly 
because of our analysis of the two leading reactions ddn and ddp. It is interesting 
in fact to study the contribution of each rate to the total error cr2H2H. The results 
are reported in Table El The role of ddn and ddp is now of the order of 50%, 
while in our previous study it was largely dominant [201 • The reduced estimated 
uncertainties on these two processes, which are very efficient burning channels, 
is mainly responsible for that. Notice that uncertainty is now dominated by a^^^. 
Even more than in the past, the role of as a baryometer is clearly established. 

- 3He- 

The estimated value for ^He almost saturates the bound reported in Section 4.1.2, 
XsHe/Xp < (l.l±0.2)-10-5^ while is a factor two smaller than the conservative 
limit, Xsiie/Xp < (1.9±0.6)-10-5. It would be of great impact for standard BBN, 
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rate 


'^^3He3He/^3He3He(%) 


heSdp 


80.7 


dp'j 


16.8 


ddp 


1.3 


ddn 


1.2 



Table 6. The contribution of reaction rate errors to the total (squared) uncertainty 
on XsHo/^p 



rate 




wnp 


98.5 


ddn 


1 


ddp 


0.25 


dn'-f 


0.25 



Table 7. The contribution of reaction rate errors to the total (squared) uncertainty 
on Yp 

as well as for theories on galactic chemical evolution, to have further experimental 
information on the abundance of this nuclide, trying to reduce the present bound. 
The error budget study shows that the theoretical uncertainty is equally shared by 
the nuclear rate errors and the uncertainty on cu^. Actually it is worth stressing 
that the value of the error due to nuclear rates is very close to what is found for 
^H. It receives the main contributions from the processes he3dp and dp^, as can be 
read from Table IHl On the other hand the fact that '^He is less critically depending 
on Ub with respect to shows up in a relatively smaller value of a^^^. 

- ^He- 

The value of the Helium mass fraction parameter is in very good agreement with 
previous theoretical determinations (see for example [1^], |2D1, |195j ). Comparison 
with the experimental values only shows a satisfactory agreement if we adopt the 
conservative estimate discussed in Section 4.1.3. The value given by jl74j (see 
Equation 4.2) is instead compatible at the 3-0" level, even considering their estimate 
of systematic error in the measurement. Similarly the result of |175j reported in 
Equation 4.3 is significantly lower with respect to the theoretical value. If we assume 
that Deuterium results points towards the consistency of both standard BBN and 
the WMAP result on ut,, we are led to the conclusion that the uncertainty of present 
measurements of Yp is largely dominated by systematic effects, which lead to an 
underestimation of ''He density. Of course a better agreement between data and 
theoretical estimates is obtained in less conservative theoretical frameworks with 
exotic physics like the presence of extra particles or large chemical potentials. These 
scenarios have been recently studied in jl94j , |195j and [20] • We show in Table [7| 
the main contribution of rate uncertainties to the squared error cr4Hc4He- We see in 
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rate 


Sal 7, /al 7, (%) 


helnxy. 


40.9 


he^heS'y 


25.1 


heldpoi 


16.2 


he3dp 


8.6 


dp-) 


4.0 


helnp 


2.7 


ddn 


1.8 


others 


0.7 



Table 8. The contribution of reaction rate errors to the total (squared) uncertainty 
on Xr-^i/Xp 



this case again how the new estimate for ddn and ddp uncertainties has lowered 
their role in the theoretical uncertainties on light nuclides. For ^He this uncertainty 
is now almost entirely due to neutron lifetime error. This is also the finding of jl9j . 

- 6Li- 

The best estimate of ^Li, Xeu/Xp = 1.1- 10^^^ gives for the ratio ^Li/'^Li the 
result 0.22-10"^, well below the result obtained for this ratio in halo stars, which, 
as discussed in Section 4.1.5, is three order of magnitude larger. Theoretical 
uncertainty on this nuclide is still very large, at the level of 100%, and is entirely 
due to the effect of ad'j, whose rate is still quite poorly known (see Section 3.4). 

- ^Li- 

There are several reactions contributing to the total uncertainty for this nuclide, 
as shown in Table |H1 Leading contributions to the 8% error on the theoretical 
estimate for uJb = 0.023 comes from be7na, he^heS'-y and be7dpa. ^Li still has 
the largest theoretical error among the observed nuclides, of the same order of 
the one induced by uncertainty. In view of this it is worth stressing that new 
experimental results in particular on the helna and he^heS'-y would be extremely 
important. However, unless unexpected new data will change the present picture of 
the main production and destruction channels for this nuclide, it seems quite hard 
to reconcile the theoretical result with the experimental data, which are smaller 
by a factor 2^3. This difference is at the level of 4cr, more severe than in our 
previous analysis in [20] because of the different treatment of nuclear rate errors, 
and is possibly an evidence for POPII star atmospheric Lithium depletion. 

4-2. 4- Results for ujh from and ^He . 

We now study how the value of Ub is constrained by standard BBN alone, without 
using the WMAP results for this parameter, but rather determining its best value and 
uncertainty from light nuclide abundances only. The aim of this analysis is to compare 
the present capability of BBN to fix the baryon content of the observable Universe with 
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respect to that of CMB anisotropies. We consider two different analysis. We first use 
the abundance only, to check the role of deuterium as a good baryometer. To this 
aim we construct the likelihood function 

, 2 



C oc exp 



(4.22) 



where notation has been described in Section 4.2.1. Using the QSO averaged result 
reported in ()4.1|) the best fit value and 68% C.L. uncertainty are found to be Ub = 
0.021±0.002. Notice that W2h2h(i^6) is dominated by the experimental uncertainty 
in QSO measurements. Improved experimental measurements of abundance would 
correspond to a more refined determination of Ub, since the role of theoretical uncertainty 
due to nuclear rate errors is presently quite small. 



1 .4 



1,2 



0,6 



0,4 



0,2 









WMAP 
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He4 


1 , 1 







0,005 0.01 0.015 0,02 0,025 0,03 0,035 0,04 



Figure 23. The uJb determination using different cosmological observables. We 
report the likeUhood function for WMAP result on CMB anisotropies (full line), 
measurement (dashed line) and ^He analysis (dotted line). 

A similar analysis can be performed using ^He. In this case we get tOb = 0.016lo;oo9, 
at 68 % C.L., which is not very constraining in view of the generous uncertainty we 
have adopted for the experimental error on Yp. 

A comparison of all results is shown in Figure ()23|) . where we report the likelihood 
profiles for the two determination of uJb from and ^He, as well as the one by WMAP. 
We read from this figure that indeed CMB anisotropies presently measure uJb with the 
best accuracy, though is also able to fix its values with a comparable precision. We 
are confident that the big effort performed by many groups in refining the theoretical 
modeling of BEN, along with new data and a better understanding of systematics may 
lead to a more accurate estimate of this important cosmological parameter by BEN 
alone. In this case comparison with present and very accurate determination from 
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CMB experiments such as PLANCK will represent a precision test of the cosmological 
model describing the evolution of the Universe. 

5. Conclusions 

In this paper we have studied in details the results of standard primordial nucleosynthesis 
in view of the fact that baryon density is now independently measured using the Cosmic 
Microwave Background anisotropics. The WMAP Collaboration result for cj^, with its 
4% accuracy, represents indeed an important step toward a detailed understanding of 
nuclide formation in the early Universe. Fixing the value of uji, in fact implies that 
standard BEN scenario becomes a parameter free model, and comparison of theoretical 
expectations with data is therefore particularly significant. 

Much effort have been done in the recent years by several groups in order to increase 
the level of accuracy of theoretical prediction on nuclide abundances, in particular by 
improving the estimate of the neutron to proton weak conversion rates, and the nuclear 
rate network. We think that the present paper represents a further contribution along 
this line, and we here summarize the main results of our study. 

We have considered in great details the n — p weak rates, which are responsible 
for neutron/proton chemical equilibrium, including the effect of distortion of neutrino 
distribution and QED effects. These rates are now known with a great accuracy, at the 
level of 0.1%. This result is now quite well established, and compatible results have 
been obtained by several groups, as we mentioned in Section 3. The main benefit of 
all these studies is that we presently have a remarkably accurate estimate of '^He mass 
fraction, which is mainly depending on the value of the neutron density at freeze out of 
weak reactions. 

We have critically reviewed the whole set of nuclear reaction rates, and we have 
reported a detailed analysis of the main rates contributing to nuclide formation, as well 
as some sub-leading channels which, however, in view of their large uncertainty, may still 
play a role in BEN. Similar studies have appeared in the recent literature [131^ which 
as the present paper largely benefit of the NACRE Collaboration database. Our main 
goal has been to get a consistent data regression method, which may be also applied to 
all cases where data shows a clear evidence for systematic errors, and which can also 
account for correlations among data belonging to the same data set without introducing 
the bias discussed in 

Our study of the main nuclear reaction rates and their corresponding uncertainties 
has been used to quantify their contribution to the total error on the estimate of light 
nuclides, in particular ^H, ^He, ^He, ^Li and ^Li. With our present regression method 
typically the role of nuclear rate uncertainty is now smaller than in previous analysis 
|14[ IT^ 1201 and is comparable or smaller than the corresponding uncertainty due to the 
present Ub error. It is remarkable that indeed for this improvement is particularly 
relevant. Refined experimental measurements of the ratio X2ii/Xp as well as of the 
baryon fraction from CME anisotropics would be a clean test of consistency of standard 
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BBN. Primordial ^Li, whose theoretical estimate is still affected by a very large error of 
the order of 100%, is not presently measured, since it is likely that the observed value 
in POPII halo stars is due to cosmic ray spallation nucleosynthesis. 

The ^He mass fraction uncertainty is estimated at the level of 0.1%, and is almost 
completely due to the neutron lifetime error. The nuclear rate uncertainties of the 
leading ddn, dip and (in7 processes contribute in fact for only 1.5% of the total 
theoretical uncertainty for uj^ = 0.023. 

The ''Li estimate is still affected by an order 10% uncertainty mainly due to the role 
of the reactions be7na, he4he3'y and heldpa. Further experimental results on these rates 
would be therefore highly desirable. In particular a better measurement of the he4he3''y 
channel would also have important implications on the estimate of solar neutrino flux. 

Apart from these processes, the BBN uncertainty budget is dominated by the effect 
of ddn, ddp and dpy, which largely contribute to both and ^He error, heSdp which 
represents the main source of uncertainty for ^He, and finally the ad'y process, which is 
entirely responsible for the large ^Li theoretical error. 

The theoretical estimates for standard BBN with experimental determinations of 
^H, ^He and ^Li show several interesting features. First of all it is remarkable that the 
value found for for the WMAP best estimate of Ub is in quite a good agreement 
with the results of QSO measurements, since they are compatible at 1 cr level. The 
value found in the same scenario for Yp is instead slightly larger than what expected 
from present data, though agreement is found at 1 a level if we adopt a conservative 
estimate of the uncertainty, as in |176j . largely dominated by systematics. The fact 
that the theoretical uncertainty for Yp is very small has important consequences in this 
respect. New experimental campaigns aimed to detect Yp and in particular a further 
effort to identify sources of systematics will clearly tell us if the standard BBN scenario 
is a satisfactory description or rather we have to consider more exotic possibilities. 

Further constraints on this may come from future experimental measurements of 
^He, which may be an independent way of fixing the baryon density. Finally, despite 
of the fact that ''Li theoretical predictions is still affected by a 10% error, nevertheless 
there is a clear indication of a possible depletion mechanism which lowers the primordial 
abundance down to the value measured with Spite plateau. A different picture might 
emerge if new experimental results on some key reaction, as for example the heldpa 
suggested in |168j . will change our present understanding of the hierarchy of the nuclear 
processes contributing to light nuclide production in the early Universe. 
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Appendix A. 

BBN reaction network 

We report in this Appendix the full list of reactions adopted in our numerical study 
of BBN responsible for the synthesis or the destruction of the A <7 nuclides. 

Weak Processes. 



Name 


Reaction 


wnp 


n p 


hS 


^ Pe + e~ + ^He 


be7 


^Be + e' ^ Ue + ^Li 



(n,7) Reactions 



Name 


Reaction 


pn'j 


p + n '^-^^ 7 + 


dn'-f 


+ n -s-^- 7 


heSn^ 


^He + n ^ 7 + ''He 


lidnj 


•^Li + n ^ 7 + ■^Li 


li7nj 


■^Li + n ^ 7 + ^Li 



(p,7) Reactions 



Name 


Reaction 


dp^ 


+ p «^ 7 + ^He 


tpi 


+ p ^ 7 + ^He 


li6p^ 


^Li + p ^ 7 + ^Be 


li7pj 


^Li + p ^ 7 + ^Be 


be7pj 


^Be + p ^ 7 + ^B 



(d,7) Reactions 



Name 


Reaction 


dd^l 


2H + 2H 4^ 7 + ^He 


adj 


^He + 2H ^ 7 + ^Li 


li7d-i 


^Li + 2h ^ 7 + ^Be 



(t,7) and (^He,7) Reactions 
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Name 


Reaction 


heSt-f 


^He + ^ 7 + ^Li 


atj 


^He + ^ 7 + ^Li 


aheS'-f 


^He + -^He ^ 7 + ^Be 


li6t'y 


6Li + '^H ^ 7 + ^Be 


helt-i 


^Be + ^ 7 + lOB 


Ii7he3-f 


^Li + =^He ^ 7 + ^°B 



{a, 7) Reactions 



Name 


Reaction 


li6a^ 


^Li + ^He ^ 7 + 


li7a^ 


^Li + ^He ^ 7 + ^^B 


be7a'y 


^Be + ^He ^ 7 + ^^C 



Charge exchange reactions 



Name 


Reaction 


heSnp 


■■^He + n ^ p + 


helnp 


^Bc + n ^ p + ''Li 



Stripping/Pickup (S/P) Reactions 



Name 


Reaction 


ddn 


+ ^ n + ^He 


ddp 


+ 2h ^ p + 


tdp 


+ 2H ^ n + ^He 


heSdp 


^He + 2H ^ p + ^He 


liddn 


^Li + 2H ^ n + ^Be 


li6dp 


6Li + 2H ^ p + ^Li 


liSpd 


^Li + p ^ 


bSnd 


8b + n ^ 2jj + 7-Qe 



and ^He S/P Reactions 
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Name 


Reaction 


li6nt 


^Li + n ^ + ^Hc 


Ii6phe3 


^Li + p ^ ^He + ^He 


heStd 


^He + ^ + ^He 


lidtd 


6Li + 3^ ^ ^ 7Lj 


Ii6he3d 


^Li + ■■^He ^ 2^ + 7gg 


liOtp 




bSnheS 


8b _^ n -s-^- ^He + ^Li 


li7tn 


^Li + ^ n + ^Be 


liVheSp 


^Li + =^He ^ p + ''Be 


heJtp 


^Be + ^ p + ^Be 


liSdt 


8Li + 2H ^ 3h 7Lj 


bSdheS 


«B + 2h ^ ^He + ^Be 


Pickup Reactions 


Name 


Reaction 


helna 


^Be + n ^ + ^He 


lilpa 


^Li + p ^ *B.e + ^He 


heOpa 


^Be + p ^ ^He + ^Li 


blOna 


lOB + n ^ + 


blOpa 


lOB + p ^ + 7Be 


lidda 


^Li + 2h ^ ^He + ^He 


h('9d(\ 


'"•Be + 2H ^He + 'Li 


belta 


^Be + ^ ^He + ^Li 


lata 


^Li + 3He ^ ^He + ^Li 


bSta 


^B + 3H ^ + 


liSheJa 


«Li + 3Hc ^ ^Hc + ^Li 


beOda 


^Be + 2h ^ ^He + ^Li 



Ingoing 2-body Reactions. 
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Name 


Reaction 


ttnn 


3h + 3h ^ n + n + ^He 


heStnp 


Sflg _|_ 3jj ^ p _|_ _|_ 4jjg 


he3he3pp 


^He + ^He p + p + "^He 


li Idncy. 




be7dpa 


■^Be + ^ p + "^He + ^He 


lilheSda 


^Li + ^He ^"^11 + ^He + "^He 


heltda 


^Be + ^ 2^ ^ 4gg ^ 4gg 


liltnna 


^Li + ^ n + n + ^He + ^He 


heltpnoi 


^Be + ^ p + n + ^He + ^He 


lilheSnpa 


^Li + 3He ^ n + p + ^He + ^He 


be7he3ppa 


^Be + ^He ^ p + p + ^He + ^He 



Ingoing 3-body Reactions. 



Name 


Reaction 


pnnd 


p + n + n^-s^n + ^H 


ppnd 


p + p + n<-i>p + 2H 


pnn'-f 


p + n + n^7 + ^H 


ppn'y 


p + p + n^7 + ^He 


pnaj 


p + n + ^He ^ 7 + ^Li 



Appendix B. 

Graphical analysis reaction rates 

A qualitative understanding of the role of each nuclear reaction to the 
creation/destruction of a definite light nuclide i can be grasped by looking at the 
temperature behavior of their contribution to the right hand side of the corresponding 
Boltzmann equation for Aj. To this end we collect in this Appendix the log-log plots of 
the combinations 

versus Tg = T/IO^A, see (12 .71) . Baryon density is chosen as ojb = 0.023. We show 
results only for those processes which contribute to change Aj at leading level or for 
some interesting sub-leading reaction. 
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Figure B2. Leading processes for production and destruction of "^H. 
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Figure B3. Leading processes for production and destruction of ^He. 




Figure B4. Leading processes for production and destruction of ^He. 
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Figure B5. Leading processes for production and destruction of ®Li. 
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Figure B6. Leading processes for production and destruction of ''Li. 
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Figure B7. Leading processes for production and destruction of ^Be. 

Appendix C. 

Fits of n ^ p weak rates 

The neutron/proton weak rates have been fitted as function oi z = me/T with 
accuracy better than 0.06% with the following expressions 

1 

u{n^p) = —e^^{-qnp/z)Y.aiz-\ 0.01 < T/MeV < 10 (C.l) 

^exp{-qpr,z)j:llibiz-' 0.1 < T/MeV < 10 
^ ^ ^ 0.01 < T/MeV < 0.1 ^ ^ 
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with 












ao = 


1 


ax = 


0.15735 


a2 = 


4.6172 


as = 


— 0.40520-10'' 


04 = 


0.13875-10'^ 


05 = 


-0.59898-10"^ 


oe = 


0.66752-10^ 


a-j = 


—0.16705-10^ 


as = 


3.8071 


Og = 


-0.39140 


aio = 


= 0.023590 


an = 


= -0.83696-10 


ai2 = 


= — 0.42095-10 


ai3 = 


= 0.17675-10 


Qnp ~ 


= 0.33979 , 




-0.62173 




0.22211-10^ 


62 = 


-0.72798-10^ 




0.11571-10^ 


64 = 


-0.11763-10^ 




0.45521-10^ 




-3.7973 




0.41266 


68 = 


-0.026210 


= 


0.87934-10-=^ 


&10 = 


-- -0.12016-10-^ 




-- 2.8602 . 



(C.3) 



(C.4) 

Appendix D. 

Fits of the S/R factors for the leading reactions 

In this Appendix we report our fits of the S or R factor for each leading reaction 
and the largest energy Em ax for which the fit is valid. 

(i) Reaction pnj : p + n ^ + '^H 

R{E) = (7.31638-10^^° + 2.35455-10"^° - 1.55683-10"^^ E 
+ 5.93351-10-^^^5 -9.25443-10-^^^2 + 6.6732-10-^^^5 

- 1.82393-10-^^ E^) cm^s-^mol-^ (D.l) 
[EMax^'^MeV] 

(ii) Reaction dp^ : + p ^ 7 + ^He 

S{E) = (0.214-10^^ + 0.556-10-^ E + 0.551-10-^ E^ 

- 0.157-10-^ MeV b (D.2) 
[EMax^'2MeV] 

(iii) Reaction ddn : "^H -\-'^H n -\- ^He 

S(E) = (0.0522 + 0.370 E - 0.196 E^ + 0.0672 E^ - 0.00885 E"^) MeV b (D.3) 
[EMax^^.lMeV] 

(iv) Reaction ddp : "^H + '^H p + 

S{E) = (0.0542 + 0.205 E - 0.0240 E^) MeV b (D.4) 
[EMax^^.lMeV] 
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(v) Reaction tdn : + ^ n + ^He 
^ 26 - 0.361 £ + 248 ^^^^ ^ 

1 I / .5-0.0479 y 
V 0.0392 J 

[EMa.^2AMeV] 



(vi) Reaction he3dp : ^He + '^H <r^p + '^He 

19.5 - 22.7 £; + 61. 8 19.5 £;3 + 4.05 

S{E) = -2 MeV b 



^ _i_ ^ £'-0.201 y 



.132 J 

[EMax^SMeV] 

(vii) Reaction he3np : ^iJe + n <^ p + '^iJ 

R{E) = (0.706-10^ - 0.149-10^° E^ + 0.521-10^° ^ - 0.239-10^^ e'^ 

+ 0.617-10" E^ - 0.449-10" E^ - 0.540-10" E^ + 0.951-10" E 
- 0.375-10" E"^) cm^s-^mol-^ 

[EMax = 1 MeV] 

(viii) Reaction be7np : ^Be + n <-> p + '^Li 

R{E) = fo.470-10^° - 0.202-10" E^ + 0.349-10" E - 0.253-10" E^ 



in o 0.109-10^° 1^1 I 
+ 0.660-10^° + cm^s-^mor^ 



317' ^ 



~ \ 0.114 

[EMax ^'2 MeV] 



(ix) Reaction aheS-f : ^He + ^He 7 + '^Be 

S{E) = 0.107-10"^ + e-°-^^2B (-0.582-10"^ - 0.606-10"^ E 

- 0.154-10-2^2) MeV b 
[EMax ^ 2.5 MeV] 

(x) Reaction lilpa : Li + p ^ ^He + '^He 

S{E) = (0.0609 + 0.173 E - 0.319 E^ + 0.217 E^) MeV b 
[EMax ^l-l MeV] 



(xi) Reaction Ii6phe3 : '^Li+p^ ^He + ^He 
S{E) = (3.158 - 1.534 £;) MeV b 
[EMax ^2.5 MeV] 
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Appendix E. 

Fits of reaction rates 

In this Appendix we collect our fits for leading and sub-leading reaction rates and 
corresponding uncertainties discussed in the paper and adopted in the numerical solution 
of BBN equation set. 

We recall that reaction rates are usually found in nuclear data catalogues in the 
form / = {av)N^''~^ where Aj is the number of ingoing nuclides (note that for Aj 7^ 2, 
(av) has only a formal meaning). This means that / has dimension {cm^'mol^^)^^~^ . 
The rates Fj (in s~^) are then obtained by multiplying the (o'f)'s by all ingoing particle 
number densities. The reverse process rates can be obtained from the forward ones 
through the following equations, which can be obtained from the detailed balance 
principle 

(i) a + b ^ c + d 

{av)ab y/J-cdJ QcQd 
where Q is the Q- value of the forward reaction. 

(ii) a + 6 -^-^ c + 7 

- I 27r J UJ ^""-'^ 

where [n-y{(Tv)c'yY'^ is the lifetime of the nucleus c for the photodisintegration 
process, and M„ is the nuclear mass unit. 

(iii) a -\- h -ir^ c -\- d -\- e 

{(yv)cde ^ / 27r X^/V M^MbM^ X^/^ ggg^ ^_q,t 

{av)ab \M^t) \MMdMe) gcgage 
In these expressions ga denotes the statistical factor 2Ja + 1, being the angular 
momentum of nuclide a. Generalization to other cases is straightfoward. These 
expressions also apply to all reactions involving identical particles, 

A„(a) + N^{h) + . . . ^ A,(c) + Na{d) + . . . , (E.4) 

apart from the replacement 

a ^- 

9a^^-^ ■ (E.5) 



A 



a- 



(i) Reaction pn^y : p + n 7 + H2 

= 44060. (^1. + 0.106597 2.75037 + 4.62949^95 

- 3.52204 Tg'^ + 1.34596 Tg^ - 0.209351 r9^) (E.6) 
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prij 



67.9521 



1.25536 - 25.8289 JTg + 618.89771, - 5043.2871, 



+ 21793.6 Tg^ - 57307.2 Tg^ + 97573.4 Tg^ - 111072. Tg2 
+ 85539.6 Tg^ - 44154.6 Tgi + 14687.1 Tg^ - 2855.93 Tg^ 



+ 247.578 Tg^ 



1/2 



(ii) Reaction dpy : H2 + p ^ 7 

1.290431 
Tg^ J 



/dP7 — ^9 



2 

"3 exp 



+ i7e3 

15.7097 + 126.82171,5 



206.50971,3 



- 721.914Tg + 2120.73T93 - 369.613 Tg^ + 173.239 Tg^ 
+ 127.838 Tgi + 100.688 Tgi - 77.3717 Tg 



exp 



0.822587 

tJ . 



(-20.6078 + 134.277 Tg5 



- 148.863^95 -651.425 71, + 1513.56 Tgt - 668.149 Tgi 
+ 690.0171,^ - 11.832 Tgi - 1.7648 Tgi - 28.933771,^) 

f 2 44793 1 / 1 

- exp <^ — '■ — — } (-10.783 + 727.882 Tg3 - 7736.03 Tg 

I T95 J ^ 
+ 32828.3 Tg - 64848.5 Tgt + 84984. Tgi - 65943. Tg^ 
+ 30784.9 Tgi - 7555.12 Tgi + 679.86 Tg^J 

(iii) Reaction ddn : '^H + n + ^He 

fddn = Tg~l exp {-Tg-i} [-1.84664-10^ + 1.22986-10^ Tgi 

- 1.3761-10^ Tgi - 6.11628-10^ Tg + 1.3329-10^ Tgi 

- 1.24333-10^ Tgi - 2.72404-10^ Tg^ + 8.52947-10^ Tgi 
+ 2.2519-10^ Tgi - 2.31204-10^ Tg^ - 294342. TgT 

+ 911550. TgT -252211. Tg^" 

Sfddn = Tg-^ exp {-Tg-i } [56996.9 - 433312. Tgi + 952341. Tg 

+ 451314. Tg - 4.57126 10^ Tgi + 5.85118 10^ Tgi 

- 1.30553 10^ Tg^ - 1.5064 10^ Tgi + 320431. Tgi 



+ 426550. Tg^ 
+ 35632. Tg^ 



16565.1 TgT - 131474. TgT 



(iv) Reaction ddp : + ^ p + ^H 



fddp = Tg~^ exp {-1.06765 Tg-i} [-5.85032-10^ + 5.23171 10^ Tgi 



1.70199-10^ To 3 



2.32242-10^ Tg - 1.18812-10^ Tg5 



+ 5.28874-10^ Tgi -9.85542-10^ Tg 
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ddp 



exp {-1.04657T9-^} (-5.74135-10*^ + 5.10713-10^ Tg 



- 1.64992-10^ To 3 



2.22795-10^ Tg - 1.11871-10^ Tgs 



+ 5.02777-10^ Tgi - 9.45483-106r92) 

- exp {-1.0657679-5} (-5.864910^ + 5.2403310^795 

- 1.7036510^795 +2.32446 10^ Tg - 1.1916710^795 
+ 5.30638 10^ Tgi - 9.88923 10^ Tg^ 

(v) Reaction tdn : + ^ n + ^He 

/tdn = ^9^5 exp {-1.34274 T9-5 - 0.233098 Tg^} (-8.11449-10^ 



2.23153-10^ Tg - 2.94397-10^ Tg^ 



1.87645-10^ 



6.05116-10^79^ + 9.51966-10^79^ - 5.29011-10^79^) 

f 0.4971161 



+ 6.22657-10^79 



-0.567854 



exp • 



Sf, 



tdn 



7g~i exp {-1.37848 7g-5 - 0.237971 7g2} (-8.70182-10^ 

+ 2.41143-10^ 79 - 3.22272-10^ 7g2 + 2.07799-10^ 7g^ 

- 6.77396-10^ 7g^ + 1.07624-10^79^ - 6.03483-10*^79^) 

+ 6.20059-10«exp|-^^:^^}79-*^-^*^°^«i 

- 7g-5 exp {-1.0508879^5 _ 0.332736792} (-4.09641-10^ 
+ 1.0649-10^ 7g - 7.15272-10^ 7g2 - 1.41552-10^ 7g^ 

+ 3.92762-10^ 79^ - 1.58174-10® 79^ + 2.12803-10^ 79*^) 

.or^.o^.^s r 0.495982- 

- 6.37982-10^ exp \ —- 



rp -0.584609 

-Lc, 



(vi) Reaction he3dp : ^He + ^ p + '^He 

fhesdp = Tg-i exp (-1.4540679"^ - 0.00623408 79^ 

+ 1.14185 

- 1.75134 

- 3.32382 

- 4.95421 

- 2.51653 



-3.13359-10' 



+ 3.10384-10® exp | 



10^795 + 1.75002-10^795 - 6.20511-10^79 
10*^795 + 5.27922-10*^795 - 1.87822-10^ 79^ 
10^795 + 2.03463-10'^ 7g5 + 6.56428-10® 7g^ 
10®7g¥ - 1.77029-10® 7g^ + 1.53089-10® 7g^ 

10^79- 

1.6191 I -0.121595 
-tg 



To 
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1 

2 



he3dp 



Tg-i exp {-1.7647379-5 - 0.0115845 Tg^} (-4.05392-10^ 



+ 1.39903 

- 2.39601 

- 7.79969 

- 8.00725 

- 4.01975 



, rr, -0.42536 



10** Tg 3 + 2.41212-10^ T93 - 6.80688-10^ Tg 
10^ Tgt +5.33313-10^ Tgt + 1.67555-10^ Tg=^ 
lO^Tgi + 3.34874-10^ Tgt + 1.33272-10*^ Tg^ 
lO^Tgf - 3.23328-10^ Tg^ + 2.52841-10^ Tg^ 

lO^TgT) 

2.75408-10^ exp 



1.78958 



- Tg-i exp {-2.04847 Tg-^ - 0.00678199 Tg^} (-4.63896-10^ 

- 2.34671-10** Tg^ + 1.6138-10^ Tgi - 3.22641-10^ Tg 

+ 5.50235-10^ Tgt - 2.96688-10^° Tgi + 8.37688-10^° Tg^ 
10^° Tgi + 2.44982-10^° Tgi +4.65934-10^° Tg^ 
10^° Tg^ + 2.36306-10^° Tg^ - 5.30274-10^ Tg^ 

lO^Tg^) 

2.75528 10** T9-°°°™^^^i exp{-i:^J^| 



- 9.56773 

- 5.17303 
+ 4.80234 



(vii) Reaction heSnp : ^He + n <-> p + 



f, 



heSnp — 



7.06494 - 4.92737 + 6.73321 Tg - 13.6597 TgS + 17.1812 Tg^ 



6.62993 Tg2 - 4.53677 Tg^ + 4.83495 Tg2 - 1.22167 Tg^ 



10*^ 



Sfhe3np = 35.1 - 50.7448 ^7^+ 31.121 Tg - 13.2721Tgi + 7.90158Tg^ 

- 3.9825 Tgi + 2.04542 Tg^ - 1.4832 Tgi + 0.511089 Tg^ 

- 0.390935 Tgi + 0.390167 Tg^ - 0.109195 Tg^ + 0.121065 Tg^ 

- 0.0907356 Tg^ + 0.0284269 Tg^ - 0.0493578 Tg^ 

„nl/2 



+ 



0.025025 Tg^l ^^^10^ 



(viii) Reaction be7np : '^Be + n <-> p + '^Li 



fi 



be7np 



6.8423-10^ - 1.4988-10^° 

3 



Tg + 1.76749-10^° Tg 



- 1.05769-10^° Tg5 + 2.6622-10^ Tg^ + 2.74476-10** TgS 

- 3.35616-10^ Tg2 + 7.64252-10^ Tgi - 5.93091-10^ Tg^ 

,7 r 0.05035181 ^_3- 
)' exp I } Tg 2 



2.28294-10^ 



Tq 



be7np 



-6.01026.10' exp ^^^^|r,-t 

Tg J 
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0.05042131 



- 2.29462-10^ expj- 



- 2.8789 10^ -3.42532 + JTq 3.09635 + JTg 



X (^8.12909 - 5.62253 ^JTg + Tg ) ( 4.24739 - 3.80751 ^JTg + Tg 
X f 1.44946 - 1.86774 jTg + Tg 



73 



(E.21) 



(ix) Reaction aheS^ : ^He + ^He ^ 7 + '^Be 
1 exp {-0.481029 Tg} 



a/ie37 



(1 + 1.1791879)' 



[0.0000461656 - 0.000460361 Tg 



0.0216009 Tg^ + 0.0696278 Tg^ + 7.34661 Tg^ 



3 , rroAaa^rr4 95 1232 Tg^ 



+ 391.131 Tg^ - 187.237Tg^ + 86.1115 Tg^ -21.6302 Tg^ 
+ 3.60069 Tg^° - 0.343228 Tg^^ + 0.0181067 Tg^^ 
- 0.000356815 Tg^^^ 



(E.22) 



5f, 



2xm 



exp 1-0.228239 Tgj , 

^ 0.0000507127- 0.000482028 Tg 

(1 + 1.31654T9)^ 

2 



0.0238316T9^ + 0.0560337T9^ + 8.40897T9^ - 106.227T9^ 
+ 434.79 Tg^ - 238.48 Tg^ + 94.7573 Tg^ - 23.7058 Tg^ + 3.80071 Tg 
- 0.370295 Tg" + 0.0199336 Tg^^ - 0.000452817Tg^^) 
exp {0.253254 Tg} 



10 



1.2857 To 



(0.0000497987 - 0.000478013 Tg 



0.0233624 Tg^ + 0.0617906 Tg'^ 



5.05896 Tg^ 



+ 418.237Tg^ - 229.349 Tg^ 



92.6387 Tq^ 



102. 196 Tg^ 
23.3706 Tg^ 



+ 3.76443 Tgi" 



0.367266 Tg" + 0.0197845 Tg^^ 



- 0.000449519 Tg^3) 



(x) Reaction li7pa : '^Li + p ^ *He + '^He 



(E.23) 



flilpa — Tg 



exp 



7.73391 



Tg- 



-8.96541-10^ + 3.86917-10^ Tg^ 



+ 4.97213-10^ Tgi - 2.58516-10^ Tg + 2.64448-10^ Tg^ 
- 1.29464-10^ Tgt - 2.68313-10^ Tg2 - 1.09411-10^ Tgi 
+ 9.98996-10^ Tgf 



(E.24) 



5fi 



lilpa 



1^ 2 

■7;Tg-^ 
2 



{6 34X73 1 
— ' ^ , I (l.64256-10^ - 1.19144-10^ T95 



+ 3.36593-10^ Tg 3 - 7.68266-10^ Tg + 1.82342-10^ Tg3 
- 1.99627-10^ Tgt + 1.24618-10^ Tg^ - 5.49787-10^ Tg5 
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+ 1.42145-10«T9t)-exp|-^:p^| (-2.99794-10^ 

+ 5.70066-10^ Tg^ + 1.33079-10^ Tgi - 7.81101-10^ Tg 
+ 1.73392-10^ Tgl - 1.82843-10^ Tgt + 1.18162-10^ Tg^ 
- 5.03632-10^ Tgi + 1.10262-10^ Tgi)' 



(xi) Reaction ad-f : ^He + ^ -f + ^Li 

fad-y = 14.82 Tg-i exp {-7.435 Tg-3} (1. + 6.572 Tg 
+ r.e-lO-^Tg^ + 2.48-10-2 Tg^) 

o^o^ 3 r 7.9041 
+ 82.8rg-2 exp I 

For Tg < 10. the lower and upper bounds result 
-.9813 + .355^7^ - .0411 Tg 



^ fad-y fad'y 



+.289 + 5.612d0e-=^^9 - 2.63 e-^^^ + .773 e-^^ 



(xii) Reaction liGpheS : ^Li+p^ ^He + ^ife 

fii6phe3^Tg--s expj-^^^^j [-7.49662-10^ + 2.05335-10^ Tg^ 

+ 3.95475-10^ Tgt - 1.94116-10^° Tg + 3.79074-10^° Tgt 
- 3.43138-10^° Tgt + 1.62629-10^° Tg^ - 3.99652-10^ Tgi 
+ 4.03339-10^ Tgt 



exp I I (4.67941-10^ - 4.18832-10® Tg 

I Tg3 J ^ 

+ 1.34221-10^^ Tgi - 1.50774-10^^ Tg - 9.45974-10^ Tgt 

+ 3.60732-10^ Tgi - 2.86104-10^^ Tg^ + 9.40736-10^ Tgi 

o 8x f 4.534591 / 7 

- 1. 15472- 10® Tg3 J -exp I (-2.15414-10^ 

- 4.59395-10® Tgt + 5.52086-10*^ Tgi - 2.14539-10^° Tg 

+ 3.82668-10^° Tgt - 3.30682-10^° Tgt + 1.51652-10^° Tg^ 

- 3.6301-10*^ Tgt + 3.58414-10® Tgt 



(xiii) Reaction at^ : '^He + ^i? <-> 7 + '^Li 



1 exp {-8.4-10-^ Tg} 
% (1 + 1.78617T9)^ 
+ 99.359 Tg2 - 989.812 Tg^ + 4368.45 Tg^ + 931.936 Tg' 



f^fy = ^ T- (0.0946142 - 4.92731 Tg 

V% (I + I.786I7T9)' ^ 
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- 391.079 Tq^ + 159.231 - 34.4076 + 3.3919 Tg 
+ 0.0175562 Tg^^ - 0.0362534 Tg^^ + 0.00311188 Tg^^ 

- 0.000087144771)^^) 



Sf, 



at-y 



2V% 



exp I -9.3-10-^X9} 

5^ 3^ (0.083877 - 4.54089X9 



(1 + 1.60171X9)'^ 
+ 96.3161 Tg"^ - 1016.55 Tg^ + 4809.48 Tg'^ - 168.102 Tg^ 
+ 208.818 Tg^ - 64.6182 Tg"^ + 10.4789 Tg^ - 0.417824 Tg^ 
- 0.0645353 Tg^° + 0.00477763 Tg^^ + 0.000200272 Tg^^ 



0.0000178642 Tc 



13 



exp {-0.0111133 Tg} 



(0.0660966 



(1 + 1.63278 Tg) 

- 3.56229 Tg + 75.1382 Tg^ - 788.241 Tg^ + 3705.89 Tg^ 

- 106.986 Tg^ + 139.556 Tg^ - 7.89845 Tg^ - 1.60357 Tg^ 

- 0.175089 Tg^ + 0.0464259 Tg^° + 0.00302332 Tg^^ 

- 0.000816826 Tg^^ + 0.0000345452 Tg^^)' 

(xiv) Reaction tpj : + p <-> 7 + ^He 

ftpj = 2.2-10^ Tg-i exp { -3.869 Tg-^} (l. + .108 TgS + 1.68dOTg 
+ 1.26 Tg + .551 Tgt + 1.06 Tgi 

^ftpj ~ 0.2 ftpj 



(xv) Reaction li7p^ : '^Li +p <-> 7 + ^Be 

fmpy = Tg--s exp { -8.62567 Tg-i - 1.13752Tg2} 



3.00142-10^ 



- 1.83661 

+ 2.02374 

+ 1.31325 

+ 2.78141 

+ 1.39934 



10^ Tg + 1.76881-10^ Tg^ - 8.47723-10^ Tg^ 

lO^^Tq'^ - 1.96501-10^" Tq''^ + 7.94528-10^ To^ 



lO^^Tg^ 

lO^Tgio 

lO^Tg^^ 



^.20935-10^ Tg^ - 9.10992-10* Tg^ 



1.07853-10^ Tg 



11 



1, 



-Tq 



exp {-5.55707T9-i - I.O4I84T92} (-25145.5 
1.58997-10^ Tg^ + 1.71826-10* Tg^ 



+ 1.07873-10^ Tg 
- 8.31031-10* Tg^ + 2.12435-10^ Tg^ - 2.87231-10^ Tg^ 



+ 2.0104-10^ Tg^ - 4.38596-10* Tg 



3.52934-10* Tq^ 



+ 2.98156-10* Tg^° - 8.89207-10^ Tg^^ + 9.98509-10^ Tg^^) 
- exp {-5.20925 Tg-5 - 1.00686 Tg^} (-14997.5 + 665017. Tg 
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- 1.08801-10^ Tg^ + 1.12999-10^ Tg^ - 5.30972-10^ Tg^ 
+ 1.32888-10^ Tg^ - 1.7653-10^ Tg^ + 1.21966-10^ Tg^ 

- 2.68716-10^ Tg^ - 2.01198-10^ Tg^ + 1.70323-10^ Tg^° 

- 5.04165-10^Tg^^ + 5.61882-10^ Tg^^)] (E.35) 

(xvi) Reaction helna : Be + n ^ ^He + ^He 



fbe7na = 2.05-10^ (1. + 3760 Tg) (E.36) 
SfbeYna ~ 0.9 Uejna (E.37) 



4.0894 



(xvii) Reaction li7dna : ''Li + ^ n + ^He + ^He 
fu-jdna = 1.71-10^X9-1 exp (-^1 + 1.49-10i°Tg-i exp 

L ig J L ig 

/ 0.0257 2.6314 4.1929 i\ 
X -2.1241 + —^ + ^ + 4.1136 T93 

V ^9 Tga Tga ) 

+ 1.66-lOi^Tg-i exp|-H:^| (E.38) 

I Tg3 J 

Sfmdna ~ 0.5 /zi7dna (E.39) 

(xviii) Reaction be7dpa : ''i?e + p + "^ife + '^He 

fbe7dpa = 1.07-10^2 ^g-i exp {-12.428 Tg"^ } (E.40) 

Sfbe7dpa ~ 0.9 fbeldpa (E.41) 

Appendix F. 

Fit of nuclide abundances 

We report a fit of tlie main nuclide abundances as functions of 77 and the number 
of effective extra degrees of freedom AA^. The fits hold for 77 corresponding to the 3 a 
range of ujb as found by WMAP, i.e. (5.48 -i- 7.12)-10~^°, while AN varies in the interval 
—3 -ir 3. The fitting function is chosen for all nuclei as follows 

\n=l 71=1 71=1 

+ j2dnX''~\^N)Aexp(j2enxA , (F.l) 

71=1 / \71.= 1 / 

where x = log^^Q (?7-10^°) and the values of the coefficients are reported in Tables IfT1|F5I 
The fit accuracy is better than 0.05% for all nuclides but for ^Li (0.2%). 
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y\ 2 U / /A ^ 


■^^^tie/ p 


p 


Xtt ; / 
' L;l/ 




1.2507-10"^ 


0.27737 


1.8737 


-1.5477-10-2 


0,2 


4.9638-10"^ 


9.7400-10-2 


0.59507 


5.2207-10-2 




-2.9574-10^^ 


-4.3242-10-3 


-0.38065 


9.6345-10-2 




-9.7061-10~3 


-2.5863-10-2 


-0.81999 


0.15855 


as 


-1.2869-10-2 


5.8552-10-2 


-0.53340 


0.21032 


Oq 


-8.9681-10-3 


0.26452 


0.56298 


0.24617 


ay 


6.7369-10-3 


0.58868 


2.3498 


0.26478 


as 


4.0429-10-2 


0.99729 


4.3744 


0.27337 



Table Fl. Coefficients of the nuclei abundance fit of Equation ljF.l|) . 





X2ll/Xp 


He/ 


yp 






1.5633-10-3 


1.1181-10-2 


7.6314-10-2 


1.7068-10-2 


b2 


7.5528-10-^ 


8.9021-10-3 


5.7666-10-2 


-6.9421-10-3 


bs 


-2.5357-10-^ 


5.6831-10-3 


2.1875-10-2 


-3.4891-10-2 


b4 


-1.2834-10-3 


2.4859-10-3 


-2.3930-10-2 


-4.2916-10-2 




-1.9679-10-3 


1.0731-10-3 


-6.1388-10-2 


-4.4713-10-2 


be 


-1.6636-10-3 


4.2430-10-3 


-5.5100-10-2 


-2.8193-10-2 


bj 


6.6829-10-^ 


1.6082-10-2 


5.4745-10-2 


1.1372-10-2 


bs 


6.6106-10-3 


4.2236-10-2 


0.36162 


7.6082-10-2 



Table F2. Coefficients bi of the nuclei abundance fit of Equation IF. Hi . 

We also report the squared error and correlations (see Equations ()4.11|) and ()4.12|) ) 
for 2h, ^He and ^Li, as function of ry and AN, valid in the same ranges for rj and AA^. 
The accuracy of these fitting expressions is better than 10 %. 

cr22H2H-10^° = 0.0151 - 0.0271 x + 0.0126 + 1.277-10-3 AA^ 



- 1.288-10-^ X AA^ + 0.81-10-^ (AA^)2 (F.2) 
^^Vne-lO' = 2.74 (F.3) 
c^7LirLi-10^° = 3-013 - 9.015 X + 6.901 x"^ + 0.116 AA^ 

- 0.190 X AA^ + 2.04-10-3 (AA^)2 (F.4) 
p2H4He = -0.108 (F.5) 
p2jj,Li = - 0.255 + 0.243 x - 0.252 x"^ + 0.018 AA^ 

- 0.015 a; AA^ - 1.1-10-^ (AN)"^ (F.6) 
p4He^Li = 0.035 (F.7) 
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X 2lT / v\ 


*-'^rle 1 p 


p 






6.3841-10"^ 


3.8084-10-^ 


2.8799-10-2 


-5.3604-10-5 


Co 


-5.2844-10"^ 


-4.2935-10-^ 


-3.6642-10-2 


2.3056-10-3 


C-x 


-1.0118-10"^ 


-6.7032-10-^ 


-5.6964-10-2 


3.2785-10-3 


C4 


-6.2060-10-5 


-3.0083-10-3 


-2.6008-10-2 


2.3728-10-3 


C5 


5.5453-10-5 


5.1127-10-3 


4.4481-10-2 


-4.3057-10-4 


C6 


1.8655-10-4 


1.2626-10-2 


0.11194 


-4.1902-10-3 


C7 


1.7096-10-4 


9.5251-10-3 


8.5888-10-2 


-6.3226-10-3 


cs 


-3.0158-10-4 


-2.1193-10-2 


-0.19410 


-1.6925-10-3 



Table F3. Coefficients Ci of the nuclei abundance fit of Equation ljF.l|) . 





X2ll/Xp 








di 


8.1121-10-^ 


6.1008-10-6 


-7.8725-10-4 


1.1108-10-4 


d2 


7.4834-10-^ 


5.4409-10-5 


1.6205-10-3 


-4.7543-10-4 


d3 


1.3879-10-5 


4.1976-10-5 


2.2747-10-3 


-5.6928-10-4 




8.4894-10-^ 


-1.3248-10-5 


9.7351-10-4 


-8.0121-10-5 


(is 


-7.5232-10-6 


-7.4104-10-5 


-1.7800-10-3 


8.7257-10-4 


de 


-2.5569-10-5 


-8.4875-10-5 


-4.2767-10-3 


1.7785-10-3 


d-j 


-2.4336-10-5 


2.4265-10-5 


-2.9394-10-3 


1.4867-10-3 


ds 


3.7481-10-5 


3.2913-10-4 


8.5262-10-3 


-2.1277-10-3 



Table F4. Coefficients di of the nuclei abundance fit of Equation (|F.1|I . 





X^-n/Xp 




yp 


Xp 


ei 


4.8889 


6.8625 


-1.8375 


1.6107 




-2.7519 


-5.9003 


5.9213 


-3.5500 


es 


-2.4691 


-5.7748 


-0.69772 


3.1784 


64 


-0.23165 


2.1014 


-4.5914 


-3.7802 




1.7868 


2.4221 


-5.3678 


5.2812 


ee 


-2.4201 


-0.75970 


5.4446 


-2.9140 



Table F5. Coefficients ei of the nuclei abundance fit of Equation (|F.1|I . 
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